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ABSTRACT 


A component of the Offiee of Naval Researeh (ONR) funded Windy Islands 
Soli ton Experiment (WISE) was eondueted from 13-15 April 2005 on the eontinental 
shelf in the northeast portion of the South China Sea to study the effects of nonlinear 
internal waves on the transmission of a 400-Hz signal. To eapture the ocean variability 
along the aeoustie path, a series of environmental moorings were deployed that sampled 
the water eolumn. Signilieant variability in the sound-speed field was observed to be 
indueed by nonlinear internal tides with a broad 10 km) horizontal scale (referred to as 
the “long-wave pattern”) and narrow (< 1 km), high-frequeney, nonlinear internal 
depression and elevation waves superimposed on the internal tides. Through the use of an 
empirieal sound-speed field and a eoupled, normal-mode aeoustie propagation model, the 
phenomenology of the nonlinear internal wave field upon the observed intensity pattern 
was examined. Analysis of the observed and modeled aeoustie intensity time-series 
indieates that the long-wave pattern dictates, to a large degree, the temporal ehanges in 
the vertieal structure of the sound intensity level. Eurthermore, both measurement and 
model results show that when the thermoeline was rapidly displaeed by the nonlinear 
internal waves, sound intensity fluctuations reached their maximum. Modeling results 
suggest that these maximums are due to the scattering of aeoustie energy into both higher 
and lower aeoustie modes along the edges of the elevation/depression waves where 
strong horizontal sound-speed gradients were present. An additional goal of this paper is 
to propose and validate an extended statistieal theory that relates the observed statistics of 
the aeoustie intensity to the number of resolvable arrivals. The number of resolvable 
arrivals depends on signal bandwidth and the eriteria of “well separateness” and was 
found to vary signifieantly as the nonlinear internal waves evolve along the transmission 
path. The theory is found to be pertinent when the temporal length of the window for 
caleulating statisties was expanded suffieiently in order to eolleet a sample population 
with the following oharaeteristies: 1) the standard deviation of the estimated number of 
arrivals is small, and 2) sufficient internal wave events are eaptured to ensure the phase 
distribution of the arrivals in the sample population is uniform. 
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I. INTRODUCTION 


A. NAVY REQUIREMENTS AND STOCHASTIC APPROACH 

Historically, the U.S. Navy has viewed sound propagation through the oeean as a 
deterministie problem. This antiquated and unrealistic approach not only hinders the 
warfighter’s eomplete understanding of the aeoustie problem, but it also impacts the 
future design of the U.S. Navy’s sonar systems. The Navy’s simplistie treatment of sound 
propagation through a medium needs replaeement to refleet a more eonventional view 
that sound propagation through the oeean should be treated stoehastieally. Although a 
stoehastie approaeh to sound propagation is not new, it differs from the deterministie 
approaeh in that it theorizes that an oeean-aeoustie wave guide is in a eontinuous state of 
flux so that a onee deterministie aeoustie signal is quiekly transformed into a fiuetuating 
signal in a relatively short distanee from its broadeast point of origin (Makris, 1996). 

B. STATISTICS OF ACOUSTIC INTENSITY IN A SHAUUOW-WATER 

ENVIRONMENT 

The Offiee of Naval Researeh (ONR) has funded several studies to investigate the 
stoehastie behavior of sound intensity in shallow-water regimes. Published physieal 
oeeanographie results from a 2001 study (Duda et al. (2004a and b)) eonducted along the 
shelf-break in the northeast portion of the South China Sea (SCS) revealed a highly 
dynamie environment where internal waves of varying spatial and temporal seales 
dominated the sound-speed field. 

Previous studies have determined that the internal wave field on the continental 
shelf is eomprised of nonlinear internal waves originating from two different sourees. 
The first souree is simply the loeal internal tides forced from the barotropie tides. In 
addition to the loeal internal tides, tidal flow over the shelf-break from the loeal internal 
tides generates loealized high-frequency, nonlinear internal waves upon the eontinental 
shelf. The seeond souree eomes from transbasin nonlinear internal tides that result from 
the barotropie tidal flow over the sills near the Bataan Islands in the Luzon Strait. 


1 



Additionally, if the flow across the sills is strong enough, as it is during “spring tides”, 
sufficient tidal energy exists within the transbasin nonlinear internal tides for the bore 
wave to transform as it propagates away from the sill into a series of narrow, large 
amplitude, nonlinear internal depression waves which Ramp et al. (2004) refer to as 
‘solitons’. Regardless of its shape, the waves associated with the nonlinear internal tides 
transit across the entire width of the South China Sea and can be observed in the internal 
wave field on the continental shelf The interaction among these wave types produced a 
complex internal wave field upon the continental shelf that would force rapid changes to 
the acoustic properties of the ocean (Ramp et ah, 2004). Additionally, nonlinear 
interactions between shoaling solitons and the bottom would force these internal waves to 
evolve from narrow, large amplitude depression waves into much broader wave forms. 
As the internal waves continue to shoal, the lower layer thins as the thermocline is driven 
toward the bottom and numerous elevation waves were observed to form behind the lead 
depression wave once the waves propagated beyond the “critical point” (i.e., where the 
thickness of the upper and lower layer is equal) (Duda et ah, 2004a and b, and Ramp et 
ah, 2004). 

In April 2005, an ONR sponsored follow-on experiment (Windy Island Soliton 
Experiment (WISE)) collected acoustic and physical oceanographic data from a series of 
moorings in a shallow-water region of the South China Sea inland of the shelf-break (see 
Section III.A for more information on the experiment). A primary focus of this field 
experiment was to measure the sound intensity fluctuations associated with a transmitted 
signal from an acoustic source and determine how various parameters within the ocean- 
acoustic wave guide influenced the statistics of measured sound intensity. 

C. DISSERTATION OBJECTIVES 

1. First Objective 

The first objective of this dissertation is to compare and contrast the sound 
intensity fluctuations resulting from a transmitted acoustic pulse through nonlinear, 
elevation internal waves against a signal transmitted through nonlinear, depression 
internal waves, in terms of phenomenology and statistics. 


2 



2 . 


Second Objective 


The second objective of this dissertation is as follows: Expanding upon previously 
established theories of the statistics of sound intensity fluctuations (Dyer, 1970 and 
Makris, 1996), propose and validate a hypothesis that relates the variance of the sound 
intensity level to the number of well-separated arrivals by incorporating critical key 
parameters such as bandwidth of the received signal (and its effect on the number of 
resolvable arrivals/degrees of freedom) as well as the temporal separation of the 
multipath arrivals. 

The outline of the paper is as follows. Section II contains a succinct synopsis of 
relevant literature, previous experiments/studies, and key concepts critical to 
understanding the results of this research. Sections III and IV correspond with the 
methodology and results from each of the two dissertation objectives. Section III begins 
with a discussion on data collection from the five environmental moorings and an 
examination of the internal wave field. With the end goal of identifying periods when 
either elevation or depression waves are present along the acoustic path as well as 
generating a temporally and spatially continuous representation of the sound-speed field, 
the mooring data is spectrally separated into a tidal and super-tidal bands followed by 
linear decomposition of the spectral bands. Using the mean sound-speed and dominant 
temporal and spatial structure of the decomposed sound-speed perturbation time-series, 
periods when elevation or depression waves occupied the transmission path were 
identified. Subsequently, the observed statistics of the acoustic intensity as these waves 
evolved along the transmission path are presented and examined. Finally, through the use 
of an empirical sound-speed field and a coupled normal-mode acoustic model. Section III 
concludes with examining the phenomenology of the observed acoustic fluctuations 
resulting from elevation or depression waves occupying the space near the source or the 
receiver. Section IV states the proposed hypothesis and how it relates to Dyer’s original 
statistical theory that related the expected statistical distribution of a noise field to the 
number of spectrally discrete noise sources. Extending this statistical theory to pulse 
propagation, the proposed theory attempts to relate the number of “resolvable arrivals” 

resulting from multipath pulse propagation to the statistics of the observed acoustic 
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intensity. Factors influencing the resolvability of individual arrivals, namely bandwidth 
of the received signal and fluctuations in the travel time, are examined and discussed. 
Four different methods of estimating the number of resolvable arrivals from the measured 
data are examined and the results are compared against the theoretical curve. The 
conclusions of this dissertation are presented in Section V. 
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II. BACKGROUND 


A. PREVIOUS STUDIES 

In late April and May 2001, an extensive experiment, the Asian Sea International 
Aeousties Experiment (ASIAEX), was conducted in the northeast portion of the South 
China Sea. Oceanographic and acoustic instruments were deployed by U.S., Taiwan and 
Singapore scientists, for the purpose of measuring acoustic fluctuations caused by ocean 
volume variability. Analysis of the environmental data revealed that the sound-speed 
fluctuations of the South China Sea shelf-slope region in April-May were dominated by 
two dynamic sources: 1) local internal tides and nonlinear waves excited from internal 
tidal flow across the nearby shelf-break, and 2) transbasin nonlinear internal tides and 
waves. 


1. Transbasin and Uocal Internal Tides 

Analysis of the ASIAEX data revealed the presence of internal tides of both 
diurnal and semidiurnal frequencies along the shelf-break in the northeastern South China 
Sea (Duda et ah, 2004a and b). Occasionally, flows induced by the local internal tides 
would lead to the formation of bores or high-frequency, nonlinear internal waves upon 
the continental shelf. 

2. Transbasin Nonlinear Internal Gravity Waves (Solitons) 

Observations along the shelf-break in the northeast portion of the South China Sea 
revealed the presence of solitons, which Duda et al. (2004a and b) describes as “large- 
amplitude, nonlinear internal gravity waves.” Although the exact generation mechanism 
is still being researched, these transbasin internal waves appear to originate from the 
Euzon Strait (Eig. 1). It has been proposed that these internal waves are excited by the 
strong tidal-current flow initiated during the spring tides across the shallow ridges of the 
southern portion of the Euzon Strait (Ramp et ah, 2004, Chiu et ah, 2004, and Duda et ah, 
2004a and b). The amplitude of these transbasin waves represent some of the largest 
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internal waves ever observed. Analysis of the ASIAEX data indieated that the solitons 
transit aeross the width of the South China Sea, moving west toward the shelf-break and 
into shallower water. 



j -6000 m 
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119 120 121 

Longitude (degrees) 


Figure 1. South China Sea overview. The boxed area near 22°N 117.5°E indicates 

where the shallow-water component of WISE (13-15 April 2005) was conducted. 
The contour lines represent isobaths at 1000-m spacing. 


As the solitons transition from deep water to shallow water, bottom influences 
and nonlinear effects begin to change their amplitude, horizontal length scale, and energy 
as they begin to shoal. Because of the shoaling effects, the amplitude of the depression 
waves decrease and their width increase as the depression waves move into shallower 
water and the extent of their amplitude becomes constrained by the bottom. The solitons 
would continue to decrease in amplitude and broaden until reaching a critical depth, at 
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which point numerous elevation waves form behind the broad, lead depression wave 
(Duda et al. 2004a and b, and Ramp et ah, 2004). Previous analysis of the ASIAEX data 
set revealed that the critieal depth eould be empirically estimated when the thiekness of 
the lower layer and the upper layer were equal. Because internal tides, as well as large 
amplitude solitons, would eontinuously foree the vertieal reloeation of the main 
thermocline, the eritieal depth, where the transformation would most likely oeeur, was 
not a fixed value but rather varied in both spaee and time. 

3. Influences of an Energetic Internal Wave Field on Acoustic 
Transmissions 

The eombination of the internal tides and the nonlinear internal waves near the 
shelf-break produee an extremely dynamie shallow-water environment. The varying 
temporal and spatial seale proeesses in a downward refraeting medium (foreing numerous 
bottom interaetions), along with the presenee of both high-frequeney nonlinear elevation 
and depression waves whieh dominate the eontinental shelf region, provided a unique 
opportunity to study the effeets of a rapidly ehanging oeean-aeoustie wave guide (sound- 
speed field) upon the statistics of the measured acoustic intensity. 

While these large amplitude internal waves were propagating through the 
ASIAEX area, moored and towed aeoustie transmitters broadeast both phase-modulated 
and frequeney-modulated aeoustie signals, spanning from 50-Hz to 600-Hz, to a 
hydrophone array positioned near the 120-m isobath. Data analysis revealed how these 
large amplitude depression waves affeeted aeoustie propagation near the shelf-break. Eor 
example, Chiu et al. (2004) demonstrated how depression internal waves bring about a 
rapid redistribution of aeoustie energy through severe mode eoupling. Orr et al. (2004) 
deseribed how rapid fluetuations in horizontal eoherenee ean lead to large fluctuations in 
array gain. And finally, Duda et al. (2004a and b) correlated scintillation index maxima 
with large numbers of high amplitude internal gravity waves. 

Although the ASIAEX data set and subsequent analysis yielded new insight on 
the effeets of large amplitude internal waves upon aeoustie propagation, the data and their 
relevanee to the Navy were limited with respeet to shallow-water aeoustie propagation. 
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The ASIAEX acoustic experiment was essentially a slope-propagation experiment as the 
relevant data were collected orthogonal to the shelf-break and limited to depths between 
120-m and 350-m. Unfortunately, insufficient acoustic or physical oceanographic data 
were collected during ASIAEX that may have provided insight into the evolution of these 
large amplitude nonlinear internal depression waves into a broad wave packet comprised 
of both depression and elevation waves as they continued to shoal into shallower water or 
how they may have affected acoustic transmissions. In April 2005, U.S. and Taiwan 
scientists attempted to rectify this “data gap” with the advent of the shallow-water 
component of the Windy Island Soliton Experiment. 

B. ARRIVAL STRUCTURES 

In the shallow-water component of WISE, the Naval Postgraduate School (NPS) 
sound source transmitted a phase-modulated signal at a carrier frequency of 400-Hz with 
a bandwidth of 100-Hz. The phase modulation employed a 5.11 second, 511 digit 
pseudo-random sequence corresponding to a compressed pulse of 10-ms resolution after 
matched-filtering. The received signals, after matched-filter processing, resulted in the 
multipath and dispersed arrivals of the compressed pulse, i.e., the pulse response (Pig. 2), 
measured as a function of depth and time (Chiu et al, 2004). The pulse response is also 
referred to as an arrival structure. Since each arrival structure represents the multipath 
spread of a single transmitted phase-modulated signal, it follows that a series of “x” 
transmissions of a phase-modulated signal will result in a sample population of “x” 
arrival structures, where “x” is simply an integer. It should be noted that an arrival may 
represents a single ray or mode, or it could be comprised of a group of ray/modes 
coherently summed. 

C. PULSE-WIDTH OF A COMPRESSED PULSE 

An inverse relationship exists between a signal’s bandwidth and the temporal 
resolution (i.e., the pulse-width) of the arrival after matched-filtering. Decreasing 
(increasing) the bandwidth of the phase-modulated signal would result in arrival with a 
wider (narrower) pulse-width (i.e., temporal resolution). Therefore, a series of arrivals 
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Pressure (Pa) 


that are broader (narrower) due to a deereased (increased) bandwidth will have a greater 
(lesser) tendency to overlap with adjacent arrivals, reducing (increasing) the number of 
temporally disjoint arrivals. 


Arrive Structure for Hydrophone 8 at 22:55:192 4/15/2CD5 (Transmission #tx527) 



Figure 2. An example of an arrival structure/pulse response collected during the 

shallow-water component of WISE. The peaks after 2.1-s represent the various 
arrivals separated in time due to multipath propagation. The arrivals may 
represent a single ray or mode, or it could be comprised of a group of ray/modes 
coherently summed. 
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D. SOUND INTENSITY (SI) AND SOUND INTENSITY UEVEU (SIU) 


For this study, the sound intensity (SI) is defined as: 


SI = - 
T 

and sound intensity level (SIL) is defined as: 


5/L = 101og 




dt 


^ T _ 

(ljuPay 


dt 


where T is the desired integration time (time differenee between the leading edge of the 
first and the trailing edge of the last arrival of interest), and p is the sound pressure of the 
signal segment in the matehed-filter output. When the output of the matched-filter is 
properly normalized, the SIL can be directly related to the transmission loss if the 
equivalent source level (SLe) is known. The equivalent source level (SLe) is defined as: 


54 -5L+ 101og “ 


\Tj 


and source level (SL) is defined as: 


5L = 101og 




{\/jPaf 


where is pulse-width of the compressed pulse (i.e., 10-ms), T is the multipath spread, 
and is the sound pressure of the compressed pulse at a range of 1-m from the source. 

With these defined values, the relationship to transform transmission loss (i.e., the model 
output) to a modeled sound intensity level is simply, SIL = SLe - TL. Note that T 
(integration time) is part of the time-bandwidth product (see next paragraph). 
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E. 


TIME-BANDWIDTH PRODUCT 


The time-bandwidth produet (TBWP) of a signal is defined as the product of the 
duration of that signal and the spectral width (bandwidth) of the medium in which the 
signal is propagating (ocean waveguide). The spectral width is simply the inverse of the 
coherence time of the ocean and the duration is the integration time used when 
calculating sound intensity (level). Coherence times of the ocean vary depending upon 
the geographic location. Typical values for the coherence time of the deep ocean, where 
low frequency, mesoscale processes dominate, are of the order of 15-30 minutes. Values 
of one-to-five minutes are more typical in shallower regions where higher frequency, 
local processes tend to determine the coherence time. 

The coherence time of the ocean along the transmission path for the shallow- 
water experiment was determined to be around 45 seconds. The integration times used to 
calculate sound intensity was less than one second. Using an assumed upper limit of one 
second for the integration time and a coherence time of 45 seconds, the spectral width is 
1/45 Hz and the TBWP was much less than unity. Following the results published by 
Makris (1996), when the measurement time (or in the case of this study, integration time) 
is much less than the coherence time, the assumption that the TBWP is equal to unity 
becomes a very good approximation. The TBWP of this study is thus taken to be unity. 

F. PRE-EXISTING STATISTICAL THEORIES ON SOUND INTENSITY 

FLUCTUATIONS 

As concluded by Makris (1996), the results of Dyer’s statistical treatment are 
limited to scenarios when the TBWP is equal to unity. The main focus of Makris’ work is 
to expand upon Dyer’s theory (1970) by considering longer integration times where the 
TBWP exceeds unity. Because the TBWP of this study does not exceed unity, the 
resulting statistical analysis of the acoustic data collected during the shallow-water 
component of WISE is expected to follow the results published by Dyer. 
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III. OBJECTIVE ONE: COMPARE AND CONTRAST THE 
PHENOMENOLOGY AND STATISTICS RESULTING EROM AN 
ELEVATION AND DEPRESSION WAVE EVENT 


In order to meet the objeetives of this study, analysis of both observed sound 
intensity measurements and modeled sound intensity were required. When portions of 
this study required the use of an aeoustic model, modeled output was compared and 
checked for consistency against in-situ measurements to ensure the model output was a 
reasonable estimate of the observed field. The specific acoustic model used in this study 
and the necessary data fields necessary to initialize the acoustic model are discussed in 
Section III.D.l. 

A. DATA COLLECTION: WINDY ISLAND SOLITON EXPERIMENT 

(WISE) 

The Windy Island Soliton Experiment was designed by Taiwan and U.S. physical 
oceanographers to observe the generation, evolution and transformation of transbasin, 
nonlinear waves in the northeastern South China Sea beginning April 2005 (Fig. 1). Until 
April 2005, no experiment had ever been specifically designed by the scientific 
community investigate the observed fluctuations resulting from transmitting an acoustic 
signal through nonlinear, elevation internal waves. WISE was delineated into two 
separate acoustic experiments, a shallow-water component (mid-April 2005), and a deep 
basin experiment (April 2005 - October 2006), respectively. 

1. Deep Basin Experiment (April 2005 - October 2006) 

Because this dissertation focuses on the data collected during the shallow-water 
component of WISE, the basin experiment will not be referred to again. However, for 
completeness, the following is provided; the objective of the deep basin experiment 
(April 2005 - October 2006) was to study and characterize the super-tidal to seasonal 
scale impacts of the transbasin, nonlinear internal waves on long range transmission loss, 
and to help monitor the evolution of the transbasin internal waves in the basin’s interior. 
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Shallow-water Component of WISE (13-15 April 2005) 


The purpose of the shallow-water component of WISE was to study the physics of 
sound propagation through nonlinear, elevation internal waves in shallow-water, and to 
quantify the associated fluctuations of the sound intensity. Beginning on April 13, 2005, 
the data collection was conducted over a three day period, prior to the deep basin 
experiment. The following is a list of the acoustic and physical oceanographic equipment 
deployed during the shallow-water component. It should be noted that analysis of the 
radar images (not shown) collected onboard the Ocean Researcher 1 indicated that the 
direction of travel for the majority of the waves in the internal wave field was roughly 
normal to the bathymetric contours (i.e., from the east-southeast toward the west- 
northwest). Thus the propagation of the acoustic signal was approximately perpendicular 
to the wave fronts of the majority of the waves in the internal wave field. 

a. NPS 400-Hz SYS13 Sound Source 

The sound source deployed for the shelf experiment was a 400-Hz source 
manufactured by Webb Research Corporation (WRC) belonging to NPS. This source has 
100-Hz bandwidth, and transmits a 5.11 second, 511 digit pseudo-random (m-sequence) 
signal at 180 dB re 1 pPa. In order to assure an adequate collection of acoustic samples, 
the sound source transmitted 30 continuous m-sequence signals lasting 153.3 seconds (or 
approximately 2 14 minutes) every five minutes during the three day deployment. 

During ASIAEX, elevation internal waves were observed to form near the 
120-m isobath. In an attempt to ensure that the transmission path would be located 
“inland” of this previously observed critical depth, the sound source was deployed along 
the 110-m isobath (Eigs. 3-4). 

b. NPS Eight Element, Vertical Line Array (VLA) 

An eight element, vertical line array (VEA) was deployed in 88-m of 
water (Eig. 4). Hydrophone spacing was eight meters between each phone resulting in an 
acoustic aperture of 56-m. With the deepest hydrophone 11.5-m above the bottom and a 

14 



surveyed bottom depth of 88-m, the depths of the hydrophones were 20.5-m, 28.5-m, 
36.5-m, 44.5-m, 52.5-m, 60.5-m, 68.5-m and 76.5-m respeetively. To faeilitate analysis 
later in this paper, the water eolumn was subdivided into three vertical sections, the upper 
portion (surface to 30-m), the middle portion (30-m to 60-m) and the lower portion (60-m 
to bottom) of the water column. The eight channel VLA digitizer sampled the acoustic 
field at 2000-Hz. 

Originally, the 400-Hz sound source and the vertical line array were to be 
separated by a distance of 20-km; however, depth restrictions of the hydrophone array 
forced the VLA to be deployed on a deeper portion of the shelf. The final distance 
between the surveyed positions of the sound source and the VLA was approximately 
17.1-km. 
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Windy Island Soliton Experiment (WISE) - Shallow Water Experiment 
April 13- 15 2005 
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Figure 3. Plane-view of the instrument orientation for the shallow-water eomponent of 
WISE. The blue, green and magenta diamonds represent the sound source 
mooring, the TidbiT temperature moorings and the VLA mooring, respectively. 
The contours represent the isobaths and are depicted every 10 meters. 
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bottom bathymetry is shown as the blue eurve. The sound souree mooring has 
been plaeed at the origin with the sound souree (blue square) positioned at a depth 
of 100 meters (10 meters off the bottom). Starting from the origin, the S, Tl, T2, 
VLA and T3 moorings have been plotted. The TidbiT sensors and the VLA 
hydrophones/TidbiT pairs are plotted as green and magenta diamonds, 
respectively. The SBE-39s are indicated by the inverted green triangle. A bottom- 
mounted Acoustic Doppler Current Profiler was deployed near the Tl mooring 
(black triangle). The horizontal axis is distance (km) and the vertical axis is the 
water depth (m). 

c. Oceanographic Data 

(1) Temperature Data. Temperature data were collected from two 
types of sensors (SBE-39 and TidbiT) located on multiple moorings. There were a total 
of four dedicated temperature moorings deployed. Starting from the vicinity of the sound 
source mooring and heading toward the VLA mooring, they were S, Tl, T2, and T3, 
respectively. Additionally, temperature sensors were placed near each hydrophone along 
the VLA mooring, providing a fifth vertical temperature profile (Eig. 4). 

The main temperature sensor used on all five temperature 


moorings were the Onset Computer Corporation’s Stow Away TidbiT temperature logger. 
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Although capable of sampling every 30 seeonds, a previous study eondueted by Woods 
Hole Oeeanographie Institution (WHOI) determined the temperature response time of the 
sensor to be approximately 10 minutes. Consequently, the TidbiT units were set to 
sample at depth at 10 minute intervals. 

More frequent temperature sampling was performed by the use of 
five SeaBird SBE-39 temperature sensors. Four SBE-39s were deployed at 60-m on the 
S, Tl, T2, and T3 mooring. The fifth SBE-39 was ahaehed to the VEA mooring and was 
deployed at 7-m (Fig. 4). The SBE-39s sampled the water column every 30 seeonds 
throughout the three day experiment. 

(2) Current Data. An upward faeing, bottom-mounted Acoustie 
Doppler Current Profiler (ADCP) was deployed next to the Tl mooring (Fig. 4). With a 
sampling rate of two samples per minute and a vertieal spaeing between samples of 4-m, 
the ADCP measured the horizontal and vertieal veloeity of the oeean eurrents during the 
three day experiment. 

B. CHARACTERISTICS OF THE OBSERVED SOUND-SPEED FIELD 

The first step in this objective was to properly identify when either elevation or 
depression waves enter and evolve along the acoustic path. Figure 5 shows the 
temperature time-series for eaeh of the five environmental moorings. The poor temporal 
resolution due to the low sampling rate of the TidbiT sensors made it diffieult to properly 
identify when a series of high frequeney elevation or depression waves were oeeupying 
the aeoustie path using the raw data. Using the teehnique published by Chiu et al. (2004), 
the temperature time-series, after being eonverted into a sound-speed time-series, was 
speetrally deeomposed into the mean sound-speed profile, and a low frequeney and high 
frequeney sound-speed perturbation time-series. The linear deeomposition of the sound- 
speed field not only aided in the proper identification of elevation/depression wave 
episodes, but the deeomposed fields were later used to develop a spatially and temporally 
eontinuous empirieal sound-speed field along the transmission path that eould be used in 
the aeoustie model. This method is discussed next. 
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Figure 5. Temperature time-series measured in Celsius (°C) reeorded by the five 

environmental moorings. The horizontal axis is time (GMT). Depth (m) is on the 
vertieal axis. 
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Empirical Orthogonal Function Decomposition of the Sound-speed 
Time-series 


Spectral densities of the sound-speed time-series calculated from the temperature 
data recorded by the SBE-39s located on the five moorings were used to aid in the proper 
selection of the two spectral bands. Spectral analysis (Fig. 6) of the three day time-series 
was expected to reveal two spectral bands, namely, tidal (i.e., internal tides) and super- 
tidal (i.e., internal waves) processes (Note that a consistent color scheme of red, blue, 
green, magenta and light blue is used in this Section to denote the data from S, Tl, T2, T3 
and the VLA moorings, respectively). The temperature time-series recorded at each of the 
five environmental mooring were converted into sound-speed time-series using an 
empirical formula (Medwin, 1975) and an approximate salinity value of 35 psu. Spectral 
peaks are clearly evident at 1, 2, 3 and 4 cycles per day (cpd) which are associated with 
the diurnal, semi-diurnal internal tides and higher order tidal harmonics observed on the 
continental shelf. The spectral energy found in the frequency bands above four cpd 
diffuses rapidly and is associated with higher frequency internal waves. Consequently, 
the tidal band was defined to have a period longer than 5 cpd, whereas the super-tidal 
band was defined to have periods shorter than 5 cpd. The use of 5 cpd to separate the two 
spectral bands is consistent with the results published by Chiu et al. (2004). It is 
interesting to note the vertical energy distribution measured by SBE-39s. Four of the 
SBE-39s were located in the lower half of the water-column (60-m) where the Brunt- 
Vaisala buoyancy frequency was observed to have a maximum. The fifth SBE-39 was 
deployed at a much shallower depth (near the top of the VLA mooring (7-m)) with the 
intention to help in the determination of mooring motion from the pressure 
measurements. As expected, the higher energy measured across all spectral bands by the 
deeper temperature sensors suggests that the sound-speed perturbation maxima for both 
bands will be found in the lower half of the water column. 
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Sound speed perturbation spectral density measured by SBE-39 



Figure 6. Spectral density of the calculated sound-speed time-series based on the data 
recorded by the five SBE-39s located on the S mooring (red), T1 mooring (blue), 
T2 mooring (green), T3 mooring (magenta) and the VLA mooring (light blue). 
Note that the SBE-39 on the VEA mooring was significantly shallower 
(positioned at 7-m) than the SBE-39s on the other four environmental moorings 
(positioned at 60-m). 
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Based on two spectral bands, the sound-speed time-series for each mooring was 
expanded into a linear summation of the mean profile, c, plus the perturbation 
components: 

c(r, z, t) = c{r, z) + Sc-j. (r, z, t) -f-(r, z,t) (1) 

where Scj and Scp^ represent the variations of the sound-speed contributed by the tidal 
and super-tidal (i.e., internal waves), respectively, in non-overlapping frequency bands. 
The mooring’s horizontal distance from the sound source is defined by r , z is depth and 
t is time. Furthermore, for each frequency band, i , and for each mooring location, r , the 
sound-speed perturbation time-series, can then be decomposed into a linear 

combination of vertical, empirical orthogonal functions (EOFs) or modes, /J'*, weighted 
by the corresponding amplitude time-series, . 

Sc. (r, Z,t) = Y,ir,t)■ (z, r) (2) 

n=\ 

The EOFs and their amplitudes were obtained using principal component analysis 
followed by projecting the data on the EOFs. Because the sampling rates differed 
between the SBE-39 sensors (sampling every 30 seconds) and the TidbiT sensors 
(sampling every 10 minutes), the data were linearly interpolated to a common sampling 
rate of 30 seconds prior to the EOF decomposition which essentially low-pass filtered the 
sound-speed time-series. Thus any internal wave with a period less then 10 minutes 
would not be accounted for in the linearly decomposed sound-speed time-series for each 
of the environmental moorings. Once the interpolated sound-speed time-series were 
generated, temporal and spatial structures of the sound-speed field were separated into 
the mean sound-speed profile and sound-speed perturbation in the tidal and internal wave 
bands using butterworth filters at each environmental mooring. 

Figure 7 shows the mean sound-speed profile for all five of the environmental 
moorings. All five sound-speed profiles are nearly identical to each other indicating a 
strong, downward refracting, range-independent mean sound-speed field along the 
transmission path. 
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(Cutoff Frequency is 5 CPD) 

S Mooring 
T1 Mooring 
T2 Mooring 
T3 Mooring 
VLA Mooring 


Sound Speed (nn/s) 

Figure 7. Mean sound-speed profde measured at the S mooring (red), T1 mooring 

(blue), T2 mooring (green), T3 mooring (magenta) and the VLA mooring (light 
blue). 
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Similar to the results found by Chiu et al. (2004), EOF deeomposition of the 
sound-speed perturbations proved to be an effeetive way to represent the vertieal 
strueture and its assoeiated temporal fluetuations within eaeh speetral band. Table 1 
shows the amount of temporal varianee eaptured by the first two modes in both the tidal 
and super-tidal bands. In the tidal band, over 92% of the temporal varianee at eaeh of the 
five environmental moorings was eaptured in just the first two modes. The higher 
frequency band was almost as good, with 73% (or higher) of the variance captured within 
the first two modes. Ignoring the VLA mooring, the principal mode (i.e., mode 1) alone 
captured over 80% (57%) of the sound-speed perturbation variance found in the tidal 
(super-tidal) band. For both the tidal and super-tidal bands, the principal mode modal 
maxima was consistently located at 60-m along the transmission path. The 60-m 
correlated well with the depth where the previously measured Brunt-Vaisala buoyancy 
frequency reaches the maximum value (Fig. 8). 

Further analysis of the EOF decomposition revealed a couple of additional 
significant findings. First, the amount of variance captured by the first mode for both 
spectral bands was range dependent. A larger portion of the temporal variance is captured 
by the principal mode from the deeper moorings (i.e., moorings S and Tl); however, 
examination of the remaining moorings reveals that the amount of temporal variance 
captured by the first mode decreases as the water depth decreases. This result suggests 
that as the waves travel toward shallower water, the bottom is modifying their structure 
and transferring energy into higher empirical modes (i.e.. Mode 2, Mode 3, etc). 
Secondly, although the first two modes from the VFA mooring and the T3 mooring 
captured 93% (73%) of the variance in the tidal (super-tidal) band, there was a distinct 
difference between the percentage of variance captured by the principal mode when 
comparing the EOF decomposition of the sound-speed perturbation field measured by the 
T3 mooring and the VFA mooring. The principal mode from the VFA mooring captured 
between 15% (2%) less than the first mode in the tidal (super-tidal) band from the T3 
mooring. This difference will be an important consideration when “tuning” the empirical 
sound-speed model. 
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Mooring 

Tidal 

Super-Tidal 

Variance Captured 

from Modes 1 & 2 


Mode 1 

Mode 2 

Mode 1 

Mode 2 

Tidal 

Super-Tidal 

S 

85.1% 

7.4% 

69.4% 

9.6% 

92.5% 

79.0% 

T1 

84.8% 

8 .8% 

68 .8% 

8.9% 

93.6% 

77.7% 

T2 

80.3% 

12.7% 

68 .8% 

9.5% 

93.0% 

78.3% 

T3 

80.0% 

12.5% 

57.3% 

15.3% 

92.5% 

72.6% 

VEA 

65.7% 

27.3% 

55.4% 

18.7% 

93.0% 

74.1% 


Table 1. Variance captured (percentage) per spectral band by the first two modes. 

In a previous study, Chiu et al. (2004) was able to use the amplitude time-series 

from his study (May 8, 2001, respectively) to isolate and track large amplitude solitons 

(i.e., individual waves) as they traveled from the deeper South China Sea basin across the 

shelf-break. Analysis of the dominant modal amplitude, with i = iw (Fig. 9, right), 

indicated that the super-tidal amplitude time-series from Chiu’s study and the super-tidal 

amplitude time-series from this study were radically dissimilar. In Chiu’s study, the large 

amplitude transbasin solitons dominated the deep water, high-frequency internal wave 

field found just off the shelf-break. Conversely, the complex, internal wave field 

observed on the shelf during the shallow-water component of WISE was composed of 

numerous smaller amplitude internal waves originating from various sources. As 

previously mentioned in Section II.A.2, it was expected (and had been previously 

observed during ASIAEX, Ramp et al. (2004)) that the observed high frequency internal 

wave field on the shelf would be a complex wave field comprised of shoaling transbasin 

solitons, the associated elevation waves forming behind the shoaling transbasin solitons, 

local internal tides and internal waves excited by internal tides flowing over the nearby 

shelf-break. All of these processes facilitate the generation of a complex, energetic 

internal wave field inside the shelf-break. Although the method of using the dominant 
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modal amplitude, af ^ time-series to traek individual elevation or depression waves along 
the aeoustie path as outlined by Chiu et al. (2004) proved unfeasible due to the eomplex 
internal wave field observed on the shelf, the EOF deeomposition was able to elearly 
identify elevation and depression wave paekets (i.e., episodes) onee the sound-speed 
time-series was reeonstrueted using equation (1) and the dominant mode for eaeh speetral 
band (Fig. 8). Analysis of this field is diseussed next. 
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EOF (Mode) 1 

(Cutoff Frequency is 5 CPD) 


Tidal Super-Tidal 




Figure 8. Principal mode (i.e., mode 1) vertical structure based on the EOF 

decomposition of the sound-speed perturbation time-series measured at each 
mooring. The principal mode for the tidal band is on the left and principal mode 
for the super-tidal band is on the right. 
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Figure 9. Spanning 13-15 April 2005, the amplitude (Principal Component) time-series 
of the principal mode associated with the tidal (left figure) and super-tidal (right 
figure) band. Starting from the sound source and continuing along the 
transmission path are the S (red), T1 (blue), T2 (green), VLA (light blue) and T3 
(magenta) moorings. 
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Identifying Periods when Elevation Waves and/or Depression Waves 
Occupy the Acoustic Path 


Examination of the mode 1 sound-speed time-series identified two periods (white 
and light blue arrow in Fig. 10, respeetively) when the thermoeline was depressed and a 
series of elevation waves were found along the transmission path. In both oases, a series 
of depression waves enter the aooustic path near the S mooring (Case 1: from 
approximately 06:00 until 15:00 on April 13; Case 2: 21:00 on 13 April until 05:00 on 14 
April). As the depression waves evolve along the acoustic path into shallower water, the 
depression waves begin to widen and compress the sound-channel toward the bottom. By 
the time the depression waves reach the T2 mooring, the depression waves have 
expanded significantly and several elevation waves have formed behind the lead 
depression wave. As the depression wave continues to move toward the VLA (or the 
nearly adjacent T3 mooring, the shoaling depression wave continues to expand 
horizontally resulting in a broad region where the sound-channel is compacted and the 
trailing elevation waves have become more pronounced. This is especially apparent when 
examining the reconstructed sound-speed time-series for the T3 mooring from 06:00 - 
13:00 on April 14 (Fig. 10). This wave transformation is consistent with the observations 
published by Ramp et al. (2004) during ASIAEX. The subsequent analysis and discussion 
will focus on the elevation waves found in Case 2 for several reasons. First, the elevation 
waves in Case 2 appear stronger in amplitude and are well defined. Second, prior to the 
entrance of broadening depression wave, a time period exists of little internal wave 
activity. The thermoeline remains relatively fixed at 40-m for several hours prior to the 
entrance of the leading edge of the depression wave observed on the S mooring at 
approximately 20:00 on April 13. This “calm period” should contribute minimal 
fluctuations to the observed sound intensity. Fastly, although the experiment was 
conducted over three days, the record length of the data collected by the environmental 
moorings did not cover the entire three day experiment due to mooring recovery. This 
abbreviated time-series placed temporal constraints on the empirical sound-speed model 
(which will be discussed in further detail in Section III.D). Case 2 occurs near the middle 
of the experiment when there is sufficient acoustic and environmental data to properly 
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examine how elevation wave packets induce fluctuations in the observed acoustic 
intensity and the phenomenology behind it. (It should be noted that the sound-speed time- 
series for the VLA mooring stops abruptly shortly after 08:00 on April 15. This “data 
gap” corresponds with the mooring recovery). 

Now that a period has been identified when elevation waves occupy the acoustic 
path, the next step is to identify a period when depression wave(s) are observed along the 
acoustic path. The depression wave observed at the S mooring around 14:00 on April 13, 
2005 was chosen for a few reasons (red arrow. Fig. 10). First, this depression wave 
maintains its structure as a depression wave as it travels along the acoustic path into 
shallower water. Unlike the transforming depression waves observed during elevation 
wave case study, this depression wave does not transform into a broader depression wave 
and elevation waves do not form behind it. Second, a calm period of minimal internal 
wave activity exists following the passage of this depression wave at each of the 
environmental moorings. As can be observed from the T3 mooring in Fig. 10, the 
depression wave exits the acoustic path at approximately 18:00 on 13 April. This occurs 
roughly the same time as the series of transforming depression waves indentified earlier 
in the elevation case study enters the acoustic path (i.e. approximately 18:00 at the S 
mooring). Finally, the depression wave occurs when there is sufficient environmental and 
acoustic data to support statistical and physical analysis of the sound intensity 
fluctuations. 

Thus an extended period has been identified (i.e., approximately from 14:00 on 
April 13 through 13:00 on April 14) where a series of elevation wave packets or 
depression wave packets occupies the transmission path. The next Section will discuss 
the statistical analysis of sound intensity fluctuations which was accomplished by 
examining the data collected by the vertical line array (VLA). The characteristics of the 
sound intensity field during this extended period will be discussed next. 
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Figure 10. Reconstructed sound-speed time-series for each of the five environmental 

moorings using only the dominant mode and its associated amplitude time-series. 
The time-series begins at 00:00, April 13, 2005 and ends at 12:00, April 15, 2005. 
The far left graph starts near the sound source (S mooring) and each consecutive 
figure follows along the transmission path (i.e., Tl, T2, VLA, and T3 mooring). 
Two periods where elevation waves form along the acoustic path have been 
identified: Case 1 (white arrow) and Case 2 (light blue arrow). Additionally, a 
depression wave transiting along the acoustic path has been identified (red arrow). 
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c. 


CHARACTERISTICS OF THE SOUND INTENSITY 


The 400-Hz sound source began its broadcast schedule at 03:00 (GMT) on April 
13, 2005. The VLA digitizer recorded the acoustic transmissions from the sound source 
for approximately the next 53 hours until the VLA mooring was recovered at 08:20 
(GMT) on April 15, 2005. Initial examination of the recorded sound pressure fields by 
the VLA digitizer revealed that the top five hydrophones experienced a sporadic 
hardware failure during the experiment. An example of this is depicted in Fig. 11. Both 
spikes and shorts could be seen in the raw voltage data. Post-experiment examination of 
the hydrophone array revealed wire damage to the array termination (which most likely 
occurred during deployment) which resulted in the intermittent failure to hydrophones 1 
thru 5. Obviously, the observed voltage spikes and shorts in the hydrophone data would 
need to be removed from the VLA data set. Otherwise, the calculated sound intensity 
measurements from the hydrophone data would contain significant errors. The procedure 
describing how these corrupted segments were removed from the raw data is discussed 
next. 

Both objectives of this dissertation focus on the characteristics of the sound 
intensity statistics. For any calculated statistical measurement, the confidence level 
associated with that value is dependent upon the number of samples used to derive it. 
Consequently, in order to ensure a high degree of confidence in the statistical analysis of 
this dissertation, only the corrupted segments were removed in order to maintain the 
largest possible sample population. Remembering that each sound source transmission 
was comprised of thirty 5.11 second m-sequences, the voltage time-series was subdivided 
into 30 segments, each with a period of 5.11 seconds. Each of these 30 segments was 
individually checked for either spikes in the voltage recordings or shorts. If a spike or 
short was found in the 5.11 second recording, the data segment was thrown out. The end 
result was removing corrupted transmission recordings while maintaining the largest 
possible sample population. With a maximum number of 30 m-sequences per 
transmission. Fig. 12 shows the percentage of good m-sequences that were recorded by 
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Amplitude Amplitude Amplitude Amplitude 


each of the hydrophones over the three day experiment. Because of the wiring damage, 
no acoustic data was collected from the majority of the middle to upper portion of the 
water column. 


Raw hydrophone data from the VLA at 04:15:19Z 04/13/2005 during transmission # tx015 




Figure 11. An example of a raw voltage plot from the eight hydrophones on the VLA. 
Spikes and shorts that occurred during one of the thirty 5.11-s m-sequence 
transmissions have been highlighted in red. These voltage segments were 
removed from the hydrophone data before sound intensity levels were calculated. 
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Figure 12. The pereentage of good m-sequenee transmissions that were reeorded on eaeh 
of the eight hydrophones. Examination reveals a signifieant laek of reliable 
acoustie data in the mid- to upper-portion of the water eolumn. 
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Figure 13. The third arrival structure as recorded by all eight hydrophones at 04:15 on 
April 13. The integration window (i.e., T) used to calculate the sound intensity 
level is outlined in red, whereas a noise segment of equal duration is outlined in 
blue. 
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Once the raw hydrophone data were reprocessed, the sound intensity levels of the 
acoustic transmissions were calculated. First, each of the 30 m-sequence segments was 
processed using a matched-filter to generate an arrival structure of the compressed pulse (i.e., 
the pulse response) (Fig. 13). The individual arrivals in the arrival stmcture correspond to the 
multipath transmissions of the compressed pulse as a function of time. An example of a 
measured arrival structure is displayed in Fig. 13. This arrival structure corresponds with the 
voltage time-series as shown in Fig. 11. Because the voltage recordings on the top three 
hydrophones were corrupted and removed, there is no associated pulse response with those 
hydrophones. The multipath spread of the signal segment (i.e., T) (bounded between the red 
lines in Fig. 13) was 0.95 seconds in length and was consistent throughout the data set. 

As previously mentioned in Section II.D, the sound intensity level (SIL) was 
defined as 


f 


5/L = 101og 


1 

r 


J p^{t)-dt 

_ 

{\pPaf 


V 


J 


where T is the multipath spread of the signal (i.e., 0.95 seconds) and p is the pressure of 
the signal segment in the matched-filter output. 

Although pulse-compressing the hydrophone data significantly increased the 
signal-to-noise ratio, the SIL estimates contain matched-filtered noise, n. Because this 
noise introduces an error in the SIL estimates, it would be useful to quantify the accuracy 
of sound intensity level estimates. In a previous study, Chiu et al. (2004) had derived the 
following equation that relates the standard deviation of the error cr^ associated with the 

SIL estimates to the standard deviation of a noise segment, cr„. 


201og(g) ^ 

V J 


In the equation above, is the sampling frequency of the signal and (T„ is the standard 
deviation of a noise segment of the same length as the multipath spread of the signal (i.e., 
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T). Using Chiu’s equation and the noise segment depieted in Fig. 13, the maximum 
standard deviation of the error in the SIL estimates for three day experiment was less than 
0.15 dB. Thus, a high level of confidence can be associated with SIL estimates and the 
resulting statistical analysis in both Objective 1 and Objective 2 of this dissertation. 

The next two Sections will discuss the analysis of both the observed sound 
intensity levels and its associated statistics during the 24-hour period identified earlier as 
a series of depression or elevation waves occupied the acoustic path. To aid in this 
discussion, the measured SIL and its associated standard deviation are presented now. 
The graph on the left in Fig. 14 is the SIL estimates time-series as measured by the 
hydrophones, while the graph on the right shows the standard deviation of the measured 
SIL time-series. Because it was desirable to examine the SIL fluctuations resulting from 
the high frequency internal wave field rather than the internal tides, the standard 
deviation of the observed SIL over a 30 minute interval (i.e., the length of the temporal 
window) was calculated. Also, as mentioned earlier, the EOF fields used to initialize the 
empirical sound-speed model (which will eventually be used to initialize the acoustic 
model) do not contain nonlinear internal waves with a period shorter than 10-mintes. Due 
to this temporal constraint in the empirical sound-speed model, the acoustic model cannot 
model sound intensity fluctuations resulting from nonlinear internal waves whose 
temporal scale is less than 10 minutes. Thus, calculating the fluctuations of the SIL using 
a 30 minute temporal window should facilitate the comparison between the observed 
variability and the modeled results. 

Finally, to prevent erroneous SIL estimates resulting from timing errors where 
only a portion of the m-sequence was contained in either the first or last 5.11-s segment 
of the hydrophone data, the first and last m-sequence for each broadcast was discarded. 
Consequently, only 28 out of the 30 possible SIL estimate “samples” contained in each 
broadcast were used when generating the statistics of this study. 

1. Vertical Structure of the SIL Time-series 

Comparison of the SIL estimates (Fig. 14, second from right) with the sound- 
speed time-series from the T3 mooring (Fig. 14, second from left) reveals that the long 
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(i.e., broad) wave pattern correlates extremely well with the SIL estimates. This suggests 
that the mean vertical distribution of the acoustic energy was driven by the longer spatial 
scale internal waves (i.e., the broad expansion and compression of the sound-channel) 
rather than the higher frequency internal waves superimposed on it. Note that when the 
sound-channel was compressed near the VLA mooring, namely, between 03:00 and 
17:00 on April 13, 02:00 and 13:00 on April 14, and finally from 03:00 until 09:00 on 
April 15, the acoustic energy was primarily confined to the lower layer of the water 
column. During these periods, the SIL estimates indicate that the energy remains focused 
near the bottom due to a compacted sound-channel near the VLA until the long-wave 
pattern expanded the sound-channel allowing acoustic energy to fill the middle portion of 
the water column once again. A good example of when the sound-channel expanded and 
acoustic energy was allowed to fill the middle portion of the water column occurred 
between approximately 19:00 on April 13 and 01:00 on April 14. During this time frame, 
the sound-channel remained expanded with the thermocline positioned at approximately 
40-m near the VLA mooring. The SIL estimates again show a broad period of time when 
the sound energy was allowed to fill the middle portions of the water column. Similarly, 
examination of the sound-speed time-series near the VLA prior to 18:00 on April 13, the 
sound-channel was compacted toward the bottom. Again, the focusing of energy near the 
bottom can be observed when examining the SIL estimate time-series prior to 16:00 on 
April 13. It is interesting to note the periodicity of the sound-channel compression near 
the VLA. Spaced roughly 24-hours apart (i.e., tidal periodicity), the sound-channel is 
routinely driven toward the bottom in the vicinity of the VLA/T3 mooring sometime 
between 02:00 and 06:00, where it remains compacted for approximately the next 12 to 
14 hours. This cyclical nature of the long-wave pattern is most likely associated with the 
internal tides. The variability of the measured SIL is discussed next. 

2. SIL Fluctuations Due to Nonlinear Elevation or Depression Waves 
near the Source or the Receiver 

Eight standard deviation maxima in the measured SIL have been identified (Fig. 
14, far right). These maxima were observed to occur when either depression or elevation 


38 



waves were in the vicinity of the sound source or the hydrophone array. The standard 
deviation maxima were observed to occur during three episodes. 

The first episode of standard deviation maxima occurs when the sound-channel is 
compacted near the receiver and a series of elevation waves transit across the hydrophone 
array (identified as events 1, 5, and 6 and on the T3 mooring sound-speed time-series 
(Fig. 14, second from the left). During all three events, the sound-channel associated with 
the long-wave pattern is compacted toward the bottom as elevation wave packets transit 
across the hydrophone array. Note that during this episode, the vertical distribution of the 
standard deviation is concentrated toward the lower half of the water column. This makes 
sense if one takes into account the compacted sound-channel associated with the long¬ 
wave pattern. As mentioned earlier, when the sound-channel is compacted, most of the 
acoustic energy will be focused toward the bottom, thus only the lowest acoustic modes 
would reach the hydrophone array as higher modes would be quickly stripped away by 
the bottom. Note, however, that the trailing and leading edge of the wave front provides a 
sharp horizontal sound-speed gradient that could allow acoustic energy trapped in the 
lowest acoustic modes to scatter into neighboring higher modes. While the elevation 
wave front may allow for the scattering of sound energy into neighboring acoustic modes, 
the compressed sound-channel would refract this energy back toward the bottom with a 
steeper incidence angle, thus quickly stripping this energy away before it was able to 
reach the hydrophones. This may explain why the measured standard deviation maxima 
of the SIL occur as the elevation waves are transiting across the hydrophone array (i.e., 
an extremely localized event). 

The second episode occurs when depression waves are near the source (identified 
as events 2, and 4 on the S mooring sound-speed time-series (Fig. 14, far left)). During 
both events, the narrow nonlinear depression waves are superimposed on the expanded 
sound-channel of the long-wave pattern. Note that the vertical distribution of the standard 
deviation is similar to that observed during the first episode in that it is concentrated 
toward the lower half of the water column. This result is expected. Unlike the previous 
episode, the expanded sound-channel associated with the long-wave pattern would allow 
for the propagation of higher acoustic modes along the transmission path. However, when 
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the depression enters the aeoustie path, the sound-ehannel is eompaeted near the souree 
thus redueing the number of aeoustie modes excited near the source. Consequently, 
higher acoustic modes would be suppressed along the acoustic path when the depression 
wave is in close proximity to the sound source. Again, with the acoustic energy limited to 
the lowest acoustic modes, any scattering (i.e., coupling) of acoustic energy along the 
acoustic path would be restricted to neighboring low modes. Consequently, because the 
acoustic energy is focused near the bottom (i.e., low modes), the greatest variability 
should occur in the lower half of the water column as acoustic energy is scattered out of 
the lowest modes. This may explain why the measured standard deviation maxima of the 
SIL are confined to the bottom of the hydrophone array. 

The final episode occurs when depression waves are in close proximity to the 
hydrophone array (identified as events 3, 7 and 8 on the T3 mooring sound-speed time- 
series (Fig. 14, second from the left). Notice that the vertical distribution of the standard 
deviation of the measured SIL is much more diffuse vertically. Again, this result is 
expected if one considers the sound-channel associated with the long-wave pattern. The 
narrow depression waves from events 3, 7 and 8 are superimposed on an expanded 
sound-channel. As mentioned earlier, the expanded sound-channel that followed these 
depression waves would allow for the excitation and propagation of higher acoustic 
modes. Unlike the previous two episodes where the acoustic energy could only be 
scattered into higher neighboring acoustic modes, the higher acoustic modes excited and 
propagating during this episode would allow for the scattering of acoustic energy into 
both higher and lower acoustic modes. Therefore, as the depression wave moves along 
the acoustic path, the trailing edge of the wave provides a sharp horizontal sound-speed 
gradient that could allow for the scattering of acoustic energy into both higher and lower 
acoustic modes, thus filling the water column with energy. This may explain why the 
measured standard deviation maxima of the SIL are more vertically diffuse during these 
events. 

It is important to summarize the results from the observed data. First, the long¬ 
wave pattern along the acoustic path associated by the expansion/compression of the 
sound-channel associated with the internal tides appears to drive the mean vertical 
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structure of the measured SIL time-series. When the long-wave pattern expands the 
sound-ehannel, the measured SIL indicates a vertically diffuse acoustic energy field 
spanning the middle and lower portions of the water column. Conversely, when the 
sound-channel is compacted toward the bottom, the acoustic energy is more focused near 
the bottom. 

Second, standard deviation maxima of the measured SIL occurred when either 
nonlinear depression or elevation waves were collocated with the sound source or the 
hydrophone array. The observed maxima in SIL fluctuations could be separated into three 
types of episodes. During the first episode, elevation wave packets were observed in the 
vicinity of the hydrophone array. With a compacted sound-channel, the standard 
deviation maxima were confined to the lower half of the water column presumably 
because the scattering of acoustic energy was limited to very low acoustic modes. 
Additionally, because the sound-channel was compacted toward the bottom, the energy 
that did scatter into higher modes was mostly likely stripped off after several interactions 
with the bottom, thus making it an extremely localized event (i.e., only when the 
elevation wave were in very close proximity to the receiver did scattered acoustic energy 
reach the hydrophone array). The events in the second episode were characterized by 
depression waves transiting over the sound source. Again, the variability in the SIL was 
confined to the lower half of the water column due to the suppression of higher acoustic 
modes near the sound source. With most of the sound energy initially limited to the 
lowest acoustic modes, acoustic energy could only be scattered out of the lowest acoustic 
modes, which consequently resulted in variance maxima observed near the lower half of 
the hydrophone array. In the final episode, standard deviation maxima occurred as 
depression waves transited across the hydrophone array. During these events, the 
expanded sound-channel associated with the internal tides presumably allowed for the 
excitation and propagation of additional acoustic modes. With a larger number of 
acoustic modes available to scatter energy, it is possible that the strong horizontal sound- 
speed gradient found along the trailing edge of these narrow, depression waves allowed 
for the scattering of acoustic energy into both higher and lower modes, thus scattering 
acoustic energy throughout the water column. 
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The next Section will discuss how a spatially and temporally continuous sound- 
speed field was generated for input into the acoustic model. This discussion will include 
what assumptions were made as well as limitations of the empirically derived sound- 
speed model. 
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EOF Decomposition of S Mooring Sound Speed 





Figure 14. The reconstructed mode 1 sound-speed time-series for the S mooring (far left) 
and the T3 mooring (second from left). Additionally, the sound intensity level 
(dB) time-series (second from right) and the standard deviation of the sound 
intensity level (dB) time-series as measured over approximately 30 minutes (far 
right). All are plotted as a function of depth (m) versus time (GMT). 
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D. MODELING RESULTS 


1. Building a Model to Understand the Phenomenology of the Sound 
Intensity Field during Nonlinear Internal Wave Events 

In order to properly interpret the phenomenology behind the observed sound 
intensity field and its assoeiated statisties, the use of an aeoustie model was required. The 
aeoustie propagation model deseribed by Chiu et al. (1995 and 1996) was used to 
examine the propagation physios of the sound intensity fiuotuations. The basio 
formulation of the model involves deoomposing the aeoustie pressure into oomplex 
envelopes that modulate (mode by mode) analytio, rapidly varying, adiabatic-mode 
solutions. Given sound-speed, density, attenuation rate, and bathymetry as a function of 
space at discrete time increments, the acoustic solution was obtained by integrating a 
coupled set of differential equations governing the coupled modal envelopes. Model 
output contained sound pressure, transmission loss, and phase and amplitude of the 
individual acoustic modes (Chiu et al., 2004). 

Chiu’s coupled normal-mode model required a sound-speed field that was 
continuous in both space (range and depth) and time. The method for generating a 
synthetic sound-speed field for ingest into the acoustic model as outlined by Chiu et al. 
(2004) was followed. In short, the empirical sound-speed model used in this study was 
constructed by advancing the principal mode’s amplitude time-series measured at the S 
mooring down-range (based on the time-step and observed wave speed in each of the 
spectral bands) at discrete time-steps and then reconstructing the sound-speed field using 
the principal mode interpolated at each range step and the corresponding amplitude value. 
Additionally, the S mooring sound-speed perturbations, , calculated with equation 2 

were slightly modified in order to account for wave decay and the dispersive nature of the 
internal waves. This was accomplished using the following relationship: 
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where fp\r,z) is the principal mode and Pp is the corresponding mode-amplitude 
time-series for spectral band i, cP is the apparent internal wave speed along the 
transmission path and a and P are the wave decay and stretching dispersion parameters, 
respectively. These are dependent on the ratio of the range-varying bathymetry, h{r) , to 
the water depth, , of the mooring used to initialize the empirical sound-speed model. 
With known bathymetry and the five environmental moorings’ principal mode, 
(r, z) , linearly interpolated along the transmission path and scaled by the local water 
depths, the model’s only free parameters were simply v and ju , and the internal wave 
speed along the transmission path for each spectral band, cP. The estimation for the 
spectral band sound-speed perturbation. Sc ., was complete once these free parameters 

were tuned to obtain a good match between the model-evolved and the observed principal 
mode amplitude time-series at each mooring. The final two parameters required by the 
model were the bottom bathymetry and the acoustic impedance of the bottom. 

A high resolution, bathymetric survey was conducted in June 2000 as part of the 
previous experiment ASIAEX. The results of this survey allowed for the extraction of a 
high resolution, three dimensional representation of bathymetry contours. A resultant 
bathymetric transect was constructed for use in the acoustic model as the lower boundary 
layer. 

Lastly, realizing the acoustic path was strongly downward refractive, reasonably 
accurate measurements of the sediment sound-speed and density were needed for model 
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initialization. A dedicated, geo-acoustic survey was not conducted during WISE; 
however, due to its close proximity to ASIAEX, bottom sediment acoustic property 
estimates from two ASIAEX studies (Marburger, 2004 and Schock, 2004) were used. 
The bottom layer was assumed to be comprised of a single homogenous layer with a 
sound-speed of 1700 m/s, a density of 1800 kg/m and an attenuation rate of 0.3 
dB/m/kHz. 

The empirical sound-speed fields generated and used in the acoustic model were 
based solely on the first mode. Because the percentage of the temporal variance captured 
by the principal mode in both spectral bands decreased concurrently with the bottom 
depth, there was concern that the primary mode (and its associated amplitude time-series) 
might well represent the sound-speed time-series measured at the deeper moorings, but 
might not fully capture the variance of the sound-speed measured at the shallower 
moorings. Additionally, because the decomposed amplitude time-series did not contain 
nonlinear internal waves with a period shorter than 10 minutes, there was concern that the 
use of solely the principal mode would not adequately capture the observed fluctuations 
in the sound-speed field. To help justify the use of the principal mode only, two figures 
were analyzed. Eirst, the sound-speed time-series using the first mode was reconstructed 
at 60-m for moorings S, Tl, T2 and T3. This time-series was then plotted against the 
measured sound-speed time-series as recorded by each of the SBE-39s (which were co¬ 
located at 60-m) (Eig. 15). The reconstructed sound-speed time-series (blue line in Eig. 
15) and the sound-speed measurements (red-dashed line in Eig. 15) show excellent 
agreement. More importantly, for all mooring ranges, the principal mode time-series 
appear to capture both lower-frequency oscillations as well as most of the higher- 
frequency internal waves from the measured data. Additionally, the data from the 
Acoustic Doppler Current Profiler (ADCP) deployed near the Tl mooring was compared 
against the vertical sound-speed profile time-series (Pigs. 16 and 17). The vertical sound- 
speed time-series for the Tl mooring has been reconstructed and plotted for April 13 and 
14. By measuring the acoustic frequency shift associated with backscatter off of 
suspended particles, the ADCP can estimate the horizontal and vertical velocities within 
the ocean. Also, by knowing the estimated sound-speed and timing between the 
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transmitted and received signal, a depth estimate can be associated with the recorded 
velocity measurements. The vertical velocity measurements from the ADCP have been 
plotted on top of the sound-speed time-series. Again, for both April 13 (Fig. 16) and 
April 14 (Fig. 17), the agreement between the in-situ velocity measurements and the 
principal mode sound-speed time-series is excellent. As depression waves moved thru 
Tl, thus forcing the thermocline toward the bottom, strong downward vertical velocity 
measurements were observed. Conversely, when elevation waves moved thru, pushing 
the thermocline toward the surface, strong upward vertical velocity measurements were 
observed. Also note the strong temporal agreement between the vertical velocity 
measurements and the sound-speed time-series, especially between 06:00 and 09:00 on 
April 14 when a series of depression and elevation waves propagated through. Based on 
the analysis of Figs. 15-17, use of only the low-passed principal mode to construct the 2- 
D empirical sound-speed field appears justified as the high frequency nonlinear internal 
waves with periods shorter than 10 minutes do not appear to contribute significantly to 
the sound-speed variability. The methodology used to determine the spectral wave speed 
and their associated tuning parameters is discussed next. 
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Figure 15. Comparison between the sound-speed time-series observed by the SBE-39s 
(red dashed line) and the EOF mode 1 reeonstructed sound-speed time-series at 
60-m (blue line) for moorings S, Tl, T2, and T3. 
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T1 Sound Speed EOF Decomposition vs In Situ Vertical Velocity for 13 April 2005 



Time (GMT) 

Figure 16. Comparison between vertical velocity measurements (black lines) as recorded 
by an ADCP co-located with the T1 mooring and the reconstructed mode 1 
sound-speed (m/s) time-series for April 13. 


T1 Sound Speed EOF Decomposition vs In Situ Vertical Velocity for 14 /^ril 2005 







Time (GMT) 

Figure 17. Comparison between vertical velocity measurements (black lines) as recorded 
by an ADCP co-located with the T1 mooring and the reconstructed mode 1 
sound-speed (m/s) time-series for April 14. 
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Re-plotting the amplitude time-series (Fig. 9) using all three-dimensions (time, 
principal component amplitude and mooring distance) made the task of tracking wave 
packets feasible (Figs. 18 and 19). Wave speed was determined by measuring the time it 
took for distinct peaks (or troughs) to travel along the acoustic path. The wave speeds 
were estimated to be 0.8 m/s and 1.22 m/s for the super-tidal and tidal bands, 
respectively. These values were consistent with values determined by Chiu et al. (2004) 
and Ramp et al. (2004). It should be noted that the wave speeds were not constant due to 
the nonlinear dispersive effect of the shoaling waves as they interacted with the bottom. 
This dispersive effect can be seen in Figs. 18 and 19. Note the peaks (or troughs) 
observed at the deepest moorings, S and Tl, are almost perfectly aligned in the vertical. 
However, as these peaks and troughs travel toward the shallower T2 and T3 moorings, 
their alignment shifts to the right. This indicates that the waves were slowing as they 
moved into shallower water. Additionally, as these waves began to shoal, the nonlinear 
bottom interactions stretched them horizontally. This was especially apparent when 
examining the peaks and troughs in the tidal band. As the waves transit from mooring S 
and finally reach mooring T3, they have been severely modified. Tuning of the empirical 
model is critical to mimic these effects. 

With the wave speed established, the only free variables left to determine were the 
tuning parameters, v and //, for each of the spectral bands. These values were 
determined through a series of iterative comparisons between the empirically propagated 
amplitude time-series originating from the S mooring and the observed amplitude time- 
series from each of the five environmental moorings (Figs. 20 and 21). The dispersion 
parameter, v, was stepped from 0 (no dispersion) to 4 (strong dispersion) at 0.1 
increments, while the wave decay parameter, //, was stepped from -3 (wave growth) to 3 
(wave decay) at 0.1 increments. Optimum values of 3.5 and 0.7 (2.7 and 1.0) were 
determined for the parameters, v and ju , for the tidal (super-tidal) band. As mentioned 
earlier, note that the empirically propagated amplitude time-series does not reach the T3 
(or the adjoining VLA) mooring until approximately 10:30 on April 13. This is the 
temporal model constraint mentioned earlier that prevented the analysis of the first 
elevation wave train observed in the principal mode sound-speed time-series (white 
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arrows, Fig. 10). Prior to 10:30, the empirically modeled internal wave field had not 
reached the hydrophones on the VLA, thus preventing its use in the acoustic model. 


Mode 1 PCs (Tidal Band) 
13-15 April 2005 



Time (GMT) 

Figure 18. Rotated view of the tidal band amplitude time-series for the S (red), T1 (blue), 
T2 (green), T3 (magenta) and VLA (light blue) moorings. Note the nonlinear 
dispersion and wave-decay as the waves evolve along the acoustic path toward 
shallower water. 
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Mode 1 PCs (Super-Tidal Band) 
13-15 April 2005 



Time (GMT) 

Figure 19. Rotated view of the super-tidal band amplitude time-series for the S (red), T1 
(blue), T2 (green), T3 (magenta) and VLA (light blue) moorings. 
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Figure 20. Comparison between the empirically propagated amplitude time-series from 
the S mooring (blue) for the tidal band and the observed principal mode amplitude 
time-series from three of the environmental moorings, Tl, T2 and T3, 
respectively (red dashed line). 
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Figure 21. Comparison between the empirically propagated amplitude time-series from 

the S mooring (blue) for the super-tidal band and the observed principal mode 
amplitude time-series from three of the environmental moorings, Tl, T2 and T3, 
respectively (red dashed line). 
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Analysis of Figs. 20 and 21 reveals that the empirieal model does not perfeetly 
mateh the in-situ data. This result was expected for a few of reasons. First, the empirical 
model used is somewhat limited in its ability to reproduce not only the complex, 
nonlinear transformation of shoaling depression waves as they broaden, but also the 
associated elevation wave packets that form behind the lead depression wave. Second, 
this model assumes that all the waves that enter the acoustic path originated from the 
same direction. This assumption facilitates the tracking of wave packets in order to 
calculate wave speed and any associated wave dispersion. Consequently, it does not take 
into account waves that enter the acoustic path from multiple directions. Analysis of radar 
images from the oceanographic vessel. Ocean Researcher 1, collected during the three 
day experiment (not shown) showed the super-tidal wave field observed on the shelf was 
also comprised of solitons reflecting off of nearby Dongsha Island. Because these 
reflected waves were propagating in a different direction than the direct path, transbasin 
waves, they would appear to transit along the environmental moorings at a different 
speed than their direct path counterparts. Thus, it would be impossible to align these 
reflected waves with the empirically propagated amplitude time-series. Finally, this 
model can only propagate wave types (i.e., elevation or depression internal waves) 
measured by the environmental mooring used for model initialization (i.e., the 
empirically propagated amplitude time-series). The empirical model can not change the 
established wave types as they propagate away from the initialization mooring. This 
limitation is examined more in the next Section. However, even with these limitations, 
there appeared to be relatively good correlation between the empirically propagated 
amplitude time-series and the amplitude time-series observed at each of the 
environmental moorings for both spectral bands. Once properly tuned, the empirical 2-D 
sound-speed fields had a temporal resolution of 30 seconds and a spatial resolution of 10- 
m and 0.5-m in the horizontal and vertical, respectively. 
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Figure 22. Three snapshots from the empirieally modeled sound-speed field along the 
aeoustie path while elevation waves were observed at the VLA. Note that 
elevation waves observed at the VLA at 09:00 whieh resulted from a compaeting 
sound-channel are not present in the modeled sound-speed field. The 400-Hz 
sound source is positioned at the origin and the VLA is approximately 17-km 
downrange. 


2. Empirical Sound-speed Model Limitations for Transforming 
Nonlinear Internal Waves 

Figure 22 shows three snapshots of the empirically modeled sound-speed field 
when elevation waves were known to transit across the hydrophone array. Beginning at 
06:00 on April 14 (top figure), the empirical sound-speed model depicts a series of 
compression and expansions of the sound-channel associated with the long-wave pattern 
along the transmission path. The broad depression of the sound-channel observed from 
approximately 3-km to 9-km downrange of the 400-Hz sound source corresponds with 
the compression of the sound-channel as measured by the S mooring from approximately 
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04:30 to 05:30 on April 14 (Fig 14). The second compression of the sound-channel in the 
empirical model was observed from approximately 14-km to 17-km downrange of the 
sound source. This compression corresponds to the broad depression of the thermocline 
observed on the S mooring at approximately 00:00 on April 14. Finally, in between the 
two sound-channel compressions, a broad area with an expanded sound-channel was 
observed in the empirical sound-speed model from approximately 9-km to 14-km 
downrange of the sound source. 

Referring back to the sound-speed time-series in Fig. 14, if the sound-speed time- 
series during an elevation wave episode was properly modeled, then the expanded sound- 
channel observed between 01:00 and 03:00 on the S mooring should compress toward the 
bottom as it propagates along the acoustic path, during which a series of elevation waves 
will form in the center of depressed sound-channel. However, comparison between the 
empirical sound-speed model and the sound-speed time-series measured by the T2 
mooring (located approximately 11.5-km downrange), indicates that by 06:00, the upper 
portion of the empirically modeled sound-channel is too shallow. The measured sound- 
speed time-series places the upper bounds of the sound-channel at approximately 60-m 
while the empirical sound-speed model positions it at approximately 40-m. Additionally, 
as observed in the measured data, while propagating from the T2 mooring to the T3A^LA 
mooring, the upper portion of the sound-channel should rapidly compress toward the 
seafloor and a series of well defined elevation waves should appear. Examination of the 
middle and lower figures in Fig. 22 clearly indicates that the expanded sound-channel 
never compressed as it propagated toward the VLA. The empirical model failed to 
reproduce the elevation wave packets observed by the environmental sensors on the 
VLA/T3 mooring. Re-examination of the empirically propagated amplitude time-series 
from Fig. 20 revealed why this may have occurred. 

Examining Eig. 20, note the two “peaks” observed in the tidal band of the 
amplitude time-series on the S mooring centered at approximately 22:00 on 13 April and 
02:00 on 14 April, which will be referred to the “early” and “late” peak, respectively. 
These two peaks correspond to two periods when the sound-channel had expanded, 
forcing the thermocline upward toward the ocean surface. Both of these periods were 
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observed on the S mooring time-series, positioned on either side of the broad depression 
wave eentered at approximately 00:00 on April 14. As diseussed earlier, the empirieal 
model propagated these peaks along the aeoustic path based on the wave speed and the 
seleeted tuning parameters. As these two peaks were propagated along the aeoustie path, 
they differed from the measured amplitude time-series observed on the downrange 
mooring in two very signifieant ways. 

First, when the “early” peak reaehed the T1 mooring (Fig. 20, seeond figure from 
the left), it eorrelates very well with the tidal band amplitude time-series from the T1 
mooring; however, when the “late” peak reaehes the T1 mooring, it does not eorrelate 
well with the measured data. The empirieally propagated “late” peak retains its 
amplitude, thus maintaining the expanding sound-ehannel, while in reality, the sound- 
ehannel had eompressed slightly between the two moorings. In-faet, eomparison between 
the model propagated amplitude time-series and the amplitude time-series measured at 
both the T2 and T3 moorings indieate that the expanded sound-ehannel assoeiated with 
the “late” peak was ineorreetly maintained and propagated along the entire length of the 
aeoustie path. 

The seeond signifieant event oeeurs when the “early” peak propagates from the 
T1 mooring to the T2 mooring (Fig. 20, middle two figures). While transiting between 
the two moorings, the “early” peak maintains its eonvex shape, while the measured data 
from the T2 mooring shows a trough at approximately 02:00 on April 14. This is 
signifieant beeause the trough observed in the T2 mooring data signifies that the sound- 
ehannel is now eompaeting toward the bottom, rather than maintaining the thermoeline 
near the surfaee (as indieated by the “early” peak). The trough is also observed to 
broaden as it propagates downrange to the T3 mooring whieh explains the broad period 
(i.e., two to three hours) when the thermoeline is depressed toward the bottom as 
observed in the measured data from the T3 mooring (Fig. 20, far right figure). The 
empirieally propagated “early” peak, however, eontinues to ineorreetly model an 
expanded sound-ehannel in the vieinity of the VLA during this time. 

This analysis showeases the empirieal model’s inability to properly model 

transforming wave types. Beeause the empirieal model was initialized using the wave 
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types from the S mooring’s amplitude time-series, it was unable to reproduce the 
evolution of the internal wave field into a series of elevation waves as observed in the 
vicinity of the VLA/T3 mooring. However, if the model was initialized using the 
amplitude time-series from the T2 mooring, the empirical sound-speed model may have 
provided a better representation of the sound-speed field as the elevation wave packets 
transformed along the acoustic path. 

With known empirical model limitation, the sound-speed fields were used within 
the acoustic model in order to generate a modeled SIL time-series. With sound pressure 
as a model output, modeled sound intensities were compared against the observed sound 
intensities and are discussed next. 

3. Examining Modeled Sound Intensity Levels during Both Elevation 
and Depression Wave Events 

With a known equivalent sound level (SLe) and range averaging the narrow-band 
transmission loss TL fields from the acoustic model into the corresponding broadband TL 
fields associated with a 100-Hz bandwidth signal (Harrison and Harrison, 1995), the 
modeled sound intensities were calculated and examined (Fig. 23). Note that an 
abbreviated time-series (i.e., approximately 24-hours) has been presented. This time 
period corresponds with the previously identified episodes when elevation or depression 
waves were collocated with either the source or the receiver. 

The modeled SIL appears to have captured several key characteristics from the 
measured SIL. First, as observed in Section III.C, the acoustic model realistically 
captured the redistribution of acoustic energy throughout the water-column associated 
with the expansion of the sound-channel at approximately 18:00 on April 13. 
Additionally, the sound-channel remained expanded until 02:00 on April 14, allowing 
acoustic energy to fill the middle and lower portion of the water column. During this 
same time-period, the modeled SIL had also expanded vertically, expanding with the 
sound-channel. Later, the sound-channel was observed to re-compress toward the bottom 
between 02:00 and 04:00 on April 14 as a series of elevation waves were observed to 
form along the acoustic path. Likewise, during this period, the acoustic model focuses the 
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sound energy toward the bottom with minimal energy eseaping into the higher modes. In 
short, the modeled SIL appears reasonably consistent with the observed intensity 
measurements. The vertical distribution of the modeled sound intensity follows the 
compression/expansion of the sound-channel which was driven by the long-wave pattern 
observed along the acoustic path. A comparison between the fluctuations observed in SIL 
estimates and the standard deviation of the modeled intensity will be discussed next. 

Figure 24 depicts the observed standard deviation of the acoustic intensity of the 
received 100-Hz bandwidth signal using a 30 minute variance window (second from the 
right) and the variance calculated from the modeled SIL (far right figure). To facilitate 
this discussion, the corresponding reconstructed mode 1 sound-speed time-series (with 
the previously identified nonlinear internal wave events) for the S mooring (figure, far 
left) and the T3 mooring (second from the left) are also presented. Finally, the periods 
previously identified when the empirical sound-speed model failed to properly replicate 
the elevation wave packets have been delineated. 

When examining the entire record length, the modeled standard deviation appears 
to correlate reasonably well with the observed variance of the received signal, with two 
exceptions. First, the variability of the modeled sound intensities when compared against 
the measured intensity fluctuations during the passage of the elevation wave packets is 
significantly weaker. As discussed earlier, examination of the empirically modeled 
sound-speed fields revealed that the elevations waves were not properly modeled. 

The second exception is that the overall magnitude of the model variance appears 
to be slightly less than the variance from the observed sound intensities. This is an 
expected result for two reasons. First, as mentioned previously, only the principal mode 
was used in the empirical sound-speed model (i.e., the temporal variability of the sound- 
speed field resulting from higher empirical modes was not captured). As such, when the 
mode 1 sound-speed fields were used with the acoustic model, the resulting modeled 
sound intensity fields would not contain the variability resulting from both the higher 
empirical modes of the sound-speed perturbation time-series as well as nonlinear internal 
waves with periods shorter than 10 minutes. Consequently, the modeled standard 

deviation time-series was expected to be less when compared against the in-situ 
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measurements. Second, because a more detailed sediment model was not available, the 
bottom was approximated as a bottom with uniform sediment properties. When 
comparing the modeled sound intensities to the measured sound intensities (Fig. 23), 
modeling the bottom as a homogeneous layer appears to be a reasonable approximation; 
however, when the real bottom was replaced with a homogeneous layer in the acoustic 
model, variability introduced to the sound intensity measurements by the real bottom was 
reduced. 

However, even with these two exceptions, the spatial and temporal variability of 
the sound intensities are reasonably well modeled. For example, from 17:00 to 21:00 on 
April 13, several standard deviation maxima in the measured SIL associated with events 
3 and 4 were observed from the hydrophone data. The model captured these events rather 
well as shown by the standard deviation maxima from 17:00 to 19:00 (Fig. 24, far right). 
Additionally, when the sound-channel was observed to collapse toward the bottom 
around 03:00 on April 14 due to the transforming/broadening depression waves resulting 
in the formation of elevation waves along the acoustic path, a standard deviation 
maximum was observed in the acoustic intensity measurements. The model captured this 
event as well with a standard deviation maximum found at approximately 03:00 as well. 

Using the results from both the empirical sound-speed model and the acoustic 
model, the final portion of this Section will conclude with a discussion relating the 
physical evolution of the measured sound intensity time-series observed during the 
passage of a depression wave along the acoustic path (i.e., events 2 & 3 identified 
earlier). 
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Figure 23. Reconstructed mode 1 sound-speed time-series for the T3 mooring (left). 

Measured sound intensity (dB) (center) versus modeled sound intensities (dB) 
(right) from 12:00 April 13, 2005 to 13:00 April 14, 2005. 


62 



































EOF Decomposition of S Mooring Sound Speed 



Figure 24. Reconstructed mode 1 sound-speed time-series for the S mooring (far left) and 
the T3 mooring (second from left). Additionally, the observed standard deviation 
of the sound intensity fluctuations (dB) (second from right) is compared against 
the modeled standard deviation (dB) (far right) using a 30-min window. The two 
periods delineate when elevation waves were observed near the receiver (which 
the empirical sound-speed model failed to replicate). 
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Examining the Phenomenology When Depression Waves Occupy the 
Acoustic Path 


Figures 25 - 27 show three modeled snapshots from the depression wave ease 
study. Eaeh of the three figures eontains the empirieally modeled sound-speed field (top), 
the modeled transmission loss XL (middle) and the relative modal magnitude from the 
aeoustie model (bottom). Figure 25 depicts the broad depression wave from events 2 and 
3 shortly after it enters the acoustic path. As previously mentioned, the downrange sound- 
channel is compacted toward the bottom between the source and the receiver. At this 
instant, the relative modal magnitude indicates that three groups of modes are initially 
excited near the sound source. The first group corresponds with the acoustic energy 
trapped by the predominately downward refracting sound-speed field found in the 
compacted sound-channel (i.e., low mode group). This mode group propagates from the 
source at low angles, with minimal interactions with the bottom. The second group of 
modes that are excited is modes 8 through 16 (i.e., intermediate mode group). These 
intermediate modes propagate away from the sound source at slightly higher angles. 
Based on their steeper trajectory and the geometry of the shallow-water experiment, the 
intermediate modes appear to interact with both the surface and bottom boundaries. 
Although the relative modal magnitude indicates that some acoustic energy from the 
intermediate modes is reaching the hydrophone array, it is significantly weaker when 
compared to the sound intensity contributed by members of the low mode group. The last 
group of modes (18-25) (i.e., high mode group) corresponds to rays that initially have 
very steep launch angles. In all three snapshots, these high angle modes are initially 
excited, but are quickly attenuated (i.e., stripped away) after several interactions with the 
bottom. 

In this first snapshot (Fig. 25), note that the sound-channel is compacted the 
length of the transmission path. With a compacted sound-channel, the majority of the 
acoustic energy was focused toward the bottom as can be observed in the transmission 
loss and relative modal magnitude figures. This agrees with the observed vertical 
distribution of the sound intensity from Fig. 23 (center). Note that at this time, only the 
lowest members (1-3) of the low mode group are excited and propagate the length of the 


64 



transmission path. Additionally, the intermediate mode group was eontributing little (if 
any) aeoustie energy to the measured SIL as it was fairly well stripped out by the bottom 
before reaching the hydrophones. The acoustic energy from the high mode groups was 
completely stripped away and did not contribute to the measured SIL. Also note that 
weak mode coupling is observed along the length of the transmission path. This first 
snapshot provides some insight into episodes 1 and 2 (i.e., elevation waves near the 
receiver and depression waves near the source, respectively). As previously inferred, the 
compression of the sound-channel has suppressed all but the lowest members of the low 
mode group. With all the acoustic energy confined to the lowest members of the low 
mode group (i.e., 1-3) the scatter of acoustic energy along the acoustic path is limited to 
neighboring higher acoustic modes. Thus, as previously concluded, the intensity 
fluctuations resulting from the scattering of acoustic energy out of the lowest members of 
the low mode group will be confined to the lower half of the water column. 

In Fig. 26, the depression wave is now centered approximately 4-km downrange 
of the sound source (top). The sound-channel remains compressed for the majority of the 
transmission path; however, as mentioned early, there is a rapid expansion of the sound- 
channel behind the depression wave. In this region between the sound source and the 
trailing edge of the nonlinear depression wave, additional members (1-6) of the low mode 
groups have been excited. Downrange of the depression wave, the sound-channel remains 
compacted, thus focusing the acoustic energy toward the bottom of the water column. 
Both the relative modal magnitude and the transmission loss plot support this conclusion. 
Note that the sharp horizontal sound-speed gradient formed during the rapid expansion of 
the sound-channel (i.e., on the trailing edge of the depression wave) is resulting in 
moderate mode-coupling, which is scattering acoustic energy into both higher and lower 
acoustic modes. The acoustic energy scattered into the lower members of the low mode 
group remains focused near the bottom due to the compacted sound-channel, while the 
acoustic energy scattered into the higher mode is stripped away approximately 11-km 
downrange, presumably due the increased number of bottom interactions associated with 
the scattering of acoustic energy into the higher modes. As observed in the previous 
snapshot, the intermediate and high mode groups are contributing little (if any) acoustic 
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energy to the SIL. The seeond snapshot provides some additional insight into the events 
eomprising the first episode (i.e., elevation waves near the hydrophone array). Previously, 
the standard deviation maxima of the measured SIL were observed to oeeur when the 
elevations were in the vieinity of the VLA. It was suspeeted that any aeoustie energy 
scattered into higher acoustic modes would be quickly stripped out as the compacted 
sound-channel strongly refracted this energy toward the bottom, thus making the 
scattering of acoustic energy into higher modes an extremely localized event. 

The final snapshot (Fig. 27) captures several of the possible processes (i.e., 
phenomenology) affecting the vertical distribution of the sound intensity time-series. 
First, note the broad expanded sound-channel behind the depression wave (top). As can 
be observed in the transmission loss figure (center) and the relative modal magnitude 
(bottom), the expansion of the sound-channel has allowed the excitation of additional 
members (i.e., 1-5) of the low mode group and acoustic energy has spread to fill the 
middle of the water column from the sound source to the trailing edge of the depression 
wave. Second, the sound-channel remains compressed for roughly the second half of the 
transmission path, thus tending to focus the acoustic energy near the bottom. However, 
strong mode-coupling along the trailing edge of the nonlinear depression wave has 
scattered acoustic energy from modes 1 thru 5 into acoustic modes 5 thru 7. Unlike the 
previous snapshot, the energy scattered into the higher acoustic modes is not stripped 
away by the bottom, but rather has filled the water column between the trailing edge of 
the depression wave and the VLA with acoustic energy. Additionally, as suspected in the 
third episode (i.e., depression waves near the receiver), strong mode-coupling is 
occurring along the trailing edge of the depression wave. Because additional members of 
the low mode group were initially excited, acoustic energy is scattered into both higher 
and lower acoustic modes, thus vertically scattering the acoustic energy throughout the 
water column. The scattering depicted in the modeling results may explain the vertically 
diffuse standard deviation maxima associated with the events comprising the third 
episode (i.e., events 3, 7 and 8). 

In conclusion, the modeling results are in good agreement with the observed 
acoustic intensity measurements and help explain its associated phenomenology. First, 
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the vertical structure of the SIL time-series follows the long-wave pattern observed along 
the transmission path. When the sound-channel was depressed (Fig. 25), the acoustic 
energy was focused toward the bottom as higher members of the low mode group were 
suppressed (i.e., episodes 1 and 2); however, when the sound-channel expanded (Fig. 27), 
the acoustic energy filled the sound-channel as additional members from the low mode 
group are excited (episode 3). Second, the rapid relocation of the thermocline (i.e., sound- 
channel compression or expansion) and its associated sharp, horizontal sound-speed 
gradient found on the trailing edge of a nonlinear depression/elevation wave was 
suspected to result in acoustic variability (episodes 1-3). Figs. 26 and 27 clearly show the 
strong mode coupling that can occur along this interface as energy is scattered into higher 
modes (episodes 1-3). And finally, even though a rapid relocation of the thermocline may 
occur, enabling acoustic energy to scatter into higher modes through mode-coupling 
along the interface, the strong, downward refracting properties that are observed during a 
compressed sound-channel may strip out the acoustic energy in a relatively short distance 
from the point of scattering, thus possibly making the associated intensity fluctuations an 
extremely localized phenomena (episode 1). 

5. Examining the Phenomenology When Elevation Waves Occupy the 
Acoustic Field near the Receiver 

As previously discussed, because it was initialized using the amplitude time-series 
from the S mooring, the empirical sound-speed model could not properly model the 
transformation of a broad depression wave into a series of trailing elevation waves; 
consequently, the acoustic model provided limited insights into the phenomenology of 
the intensity statistics as related to mode-coupling when elevation waves occupy the 
acoustic field near the receiver. However, the analysis of the acoustic modeling results in 
the previous section have provided some insight into the possible phenomenology of the 
distribution of the acoustic intensity when elevation waves pass by the hydrophone array. 
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Figure 25. Combined modeled results from the empirical sound-speed model (top), the 
modeled transmission loss (middle) and the relative modal magnitude from the 
acoustic model (bottom) as a depression wave moves along the acoustic path. 

This first snapshot depicts the center of the depression wave at approximately 2- 
km with the trailing edge of depression wave located above the sound source (i.e., 
at the origin) (top figure). 
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Modeled 2D Sound Speed Field: Tinne = Apr/13/2005 15:21:00Z 
Tidal Band: Mu = 3.5 Nu = 0.7 
Super-Tidal Band: Mu = 2.7 Nu = 1 
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Figure 26. Combined modeled results from the empirical sound-speed model (top), the 
modeled transmission loss (middle) and the relative modal magnitude from the 
acoustic model (bottom) as a depression wave moves along the acoustic path. 
This second snapshot depicts the depression wave centered approximately 4-km 
downrange from the acoustic source (top figure). 
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Figure 27. Combined modeled results from the empirical sound-speed model (top), the 
modeled transmission loss (middle) and the relative modal magnitude from the 
acoustic model (bottom) as a depression wave moves along the acoustic path. 
This third and final snapshot depicts the depression wave centered approximately 
11-km downrange from the acoustic source (top figure). 
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IV. OBJECTIVE TWO: PROPOSE AND VALIDATE AN 
EXPANDED STATISTICAL HYPOTHESIS RELATING THE 
MEASURED VARIANCE TO THE NUMBER OF RESOLVABLE 

ARRIVALS 


In Section I.C. of his paper, Dyer (1970) proposed that the statistics of an ambient 
noise field could be reasonably predicted if one considered that the acoustic field was 
generated by multiple distant, spectrally separated acoustic sources. The received 
multipath signal from each individual noise source, whose phase was uniformly random, 
would produce mean-square pressure fluctuations that would follow an exponential (i.e.. 
Chi-squared with two degrees of freedom) distribution. As the number of distant 
“sources” increases, the following will occur. First, for each additional spectrally 
separated noise source, the degrees of freedom will increase by two. Second, the statistics 
of the mean-square pressure field from the noise sources will approach a Gaussian 
distribution as the number of noise sources becomes large. And finally, as the degrees of 
freedom for the Chi-squared distribution increases, the variance of the noise fluctuations 
will decrease. In his paper. Dyer presented a set of equations (i.e., his equations 6, 28 and 
29) relating the standard deviation of the noise intensity level (dB) to the number of 
sources, specifically: 




10 p 

In(lO) VT 


S,il) 


where 


5^(1) = 0 or 5,(0 

H.=l W 

and / is defined as the number of spectrally separated noise sources. This relationship has 
been plotted in Fig. 28. 
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standard Deviation of a Phase Pandom Signal based on the Number of Disorete Souroes 



Number of Souroes 

Figure 28. Expected standard deviation for a phase-random signal based on the number 
of spectrally separated noise sources (Dyer, 1970). 

Realizing that each arrival in any given arrival structure is statistically similar to 
Dyer’s “individual distant noise source,” it is hypothesized that the density of the sound 
intensity fluctuations will follow the log of a Chi-squared distribution where the degrees 
of freedom are determined by the number of resolvable arrivals (i.e.. Fig. 2) thus allowing 
the prediction of SIL/TL variability based on the number of resolvable arrivals. It is 
further postulated that as the number of resolvable arrivals increases, the measured 
variance of the sound intensity distribution will follow the expected changes in the 
statistical distribution as predicted by Dyer’s statistical theory of a noise field (1970). The 
proposed hypothesis states that the degrees of freedom will be determined by twice the 
number of resolvable arrivals. Similar to Dyer’s theory, as the number of resolvable 
arrivals increases, the variance of the sound intensity is expected to decrease as the 
distribution approaches the log of a Gaussian distribution. The relationship between the 
measured standard deviation of the acoustic intensity and the number of resolvable 
arrivals as specified in the proposed hypothesis are expected to follow the curve shown in 
Fig. 28. Furthermore, it is hypothesized that resolvability of the arrivals is sensitive to 
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both the bandwidth of the reeeived signal and the spaeing of the multipath travel times. In 
short, eaeh arrival is hypothesized to represent an independent intensity and the number 
of resolvable arrivals (i.e., degrees of freedom) drives the statisties. The next Seetion 
diseusses two aeoustie parameters that ean affeet the resolvability of an arrival. 

A. PARAMETERS INFLUENCING ARRIVAL RESOLVABILITY 

1. Temporal Variability Induced by the Transmission Medium 

Ineessant, small-seale ehanges in the sound-speed field along the transmission 
path due to miero-seale proeesses result in a eontinuously ehanging multipath 
propagation. Beeause the transmission path is in a eonstant state of flux, the travel time 
(the time it takes for the pulse signal to propagate from an aeoustie souree to a reeeiver) 
is ehanging as well. Therefore, due to their fluetuating travel time, individual arrivals 
depleted on an arrival strueture will oseillate about a “mean travel time.” Two adjaeent 
arrivals may be well separated and show up as two distinet independent pulses on an 
arrival strueture. However, if adjaeent arrivals experienee a large enough variation in 
their travel time and are spaeed within a elose proximity to one another, then there may 
be times when the adjaeent arrivals overlap one another. This overlap will lower the 
number of resolvable arrivals, reduee the degrees of freedom and alter the statisties. 

2. Bandwidth Sensitivity 

An additional parameter that will influenee the resolvability of an individual 
arrival is its pulse-width (whieh is inversely proportional to the bandwidth of the reeeived 
signal). If the pulse-width (i.e., the temporal resolution of the pulse signal) is allowed to 
inerease, broad individual adjaeent arrivals will have a greater tendeney to overlap and 
will be more diffieult to resolve than narrower ones. This would also reduee the degrees 
of freedom and alter the statisties. Therefore, this suggests that the statisties of the sound 
intensity fluetuations, whieh is hypothesized to depend upon the number of resolvable 
arrivals (i.e., degrees of freedom), is sensitive to both the degree of the travel time 
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fluctuations and the pulse-width of the arrival. Both modeled and in-situ measurements 
were analyzed for bandwidth sensitivity and are discussed next. 

B. STATISTICAL ANALYSIS OF THE MODELED INTENSITY 

The modeled transmission loss was based on a continuous wave signal and would 
first need to be transformed into the eorresponding transmission loss fields associated 
with a broadband source. Given the fractional bandwidth of the known (or desired) 
broadband signal, and the distance separating the source and receiver, the model output 
was transformed into a broadband signal using a range averaging teehnique outlined by 
Harrison and Harrison (1995). Consequently, transmission loss fields for a modeled 10-, 
25-, 50-, 75-, 100-, 200- and 400-Hz bandwidth signal were ealeulated (not shown). Once 
the TL fields for the seven modeled broadband sourees were ealeulated, the intensity 
fluctuations from the modeled transmission loss fields were analyzed (Figs. 29 and 30). 

1. Sensitivity of the Statistics to the Bandwidth of the Modeled Pulse 

Figure 29 depicts the modeled standard deviation calculated over a 30 minute 
window for the seven modeled broadband signals. The statistical results from the 10-Hz 
bandwidth signal are presented in the upper graph followed by the results from the 
modeled 25-, 50-, 75-, 100-, 200- and 400-Hz bandwidth signal respeetively. By 
examining this figure, two important observations can be made. First, note the inverse 
relationship between bandwidth of the received signal and the magnitude of the standard 
deviation of the modeled signal. Speeifically, the magnitude of the standard deviation 
decreases as the bandwidth inereases. Even when the variance window was extended to 
180 minutes, this relationship still exists (Fig. 30). This model result is in eompete 
agreement with the proposed hypothesis. Seeond, the standard deviation values are non¬ 
stationary (i.e., time dependent). By choosing such a temporally short variance window, 
the intensity fluctuations depicted in Fig. 29 are known to be assoeiated with the high 
frequeney internal wave field as discussed in Section III. Even when the variance window 
is extended to 180 minutes, the standard deviation values continue to be time variant, 
suggesting that these values are still being influenced by the time-varying, high frequeney 
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internal wave field. Diseussed in detail in the next Section, this finding will be an 
important consideration when attempting to estimate the degrees of freedom associated 
with the measured data. 


Standard Deviation of the Modeled SIL (dB) measured over 6 Transmissions (~30 min). 
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Figure 29. The standard deviation of the modeled sound intensity (dB) associated with a 
10-, 25-, 50-, 75-, 100-, 200- and 400-Hz bandwidth signal. The standard 
deviation was calculated from the modeled broadband TL using a 30 minute 
window. 
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standard Deviation of the Modeled SIL (dB) measured over 36 Transmissions (^180 min). 
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Figure 30. The standard deviation of the modeled sound intensity (dB) associated with a 
10-, 25-, 50-, 75-, 100-, 200- and 400-Hz bandwidth signal. The standard 
deviation was calculated from the modeled broadband TL using a 180 minute 
window. 


C. STATISTICAL ANALYSIS OF THE MEASURED INTENSITY 

The 100-Flz bandwidth signal collected on the VLA was bandpass filtered in 
order to generate a 10-, 25-, 50-, and 75-Hz bandwidth signal, thus providing acoustic 
intensity measurements from five broadband sources. When the bandwidth of the original 
broadband signal was narrowed through the use of a butterworth bandpass filter, the 
pulse-width of the recorded signal was increased. This effect can be easily observed when 
the resulting arrival structures from the bandpass filtered signal are plotted (Fig. 31). The 
intensity statistics of these five broadband signals will be discussed next. 
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Sensitivity of the Statistics to the Bandwidth of the Received Pulse 


The effect of decreasing the bandwidth of the received signal to the temporal 
resolution of the arrival (i.e., pulse-width) is apparent in Fig. 31. Narrowing the 
bandwidth of the received signal will result in fewer individual arrivals. For example, 
note the two separated arrivals located at approximately 19.36-s and 19.38-s associated 
with the 50-, 75- and 100-Hz bandwidth signal. When the 100-Hz bandwidth signal was 
filtered to produce a 25-Hz bandwidth signal, the two arrivals broadened and merged into 
a single arrival. Furthermore, when the signal was ultimately filtered to a 10-Hz 
bandwidth signal, the once identifiable individual arrivals across the 25-, 50-, 75- and 
100-Hz arrival structures have been lumped together into a single broad arrival. Thus, the 
direct relationship between the number of arrivals and the bandwidth of the received 
signal is once again observed and in strong agreement with the proposed hypothesis. In 
order to facilitate a direct comparison with the statistics of the modeled acoustic 
intensities, the standard deviation of the received signals from the five broadband sources 
was calculated using a 30 minute and 180 minute variance window and is discussed next. 
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Figure 31. Starting with the 100-Hz bandwidth arrival structure (black), the broadband 
signal was bandpass filtered to derive the arrival structures for a 10-Hz (light 
blue), 25-Hz (green), 50-Hz (red), and 75-Hz (blue) bandwidth signal. 

Examination of both Fig. 32 and Fig. 33 reveals that the measured standard 
deviation of the sound intensity exhibits a bandwidth sensitivity similar to that depicted 
in the model results from the previous Section. Namely, the magnitude of the standard 
deviation decreases as the bandwidth of the broadband signal increases regardless of the 
variance window used. Also consistent with the statistical analysis of the modeled 
intensity. Figs. 32 and 33 suggest that the statistics of the measured intensity are time- 
variant and are being driven by the high frequency, internal wave field. Additionally, 
there appears to be a tidal-scale periodicity associated with the standard deviation 
maxima, suggesting that the high-frequency internal wave field is periodic. Two standard 
deviation maxima were observed 24-hours apart at approximately 18:00 on both April 13 
and April 14. It is interesting to note that during ASIAEX, transbasin solitons were 
observed to reach the 120-m bathymetric curve around 16:00 each day. Considering the 
close proximity of the shallow-water component of WISE to the ASIAEX site, the 
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standard deviation maxima observed at 18:00 on both days may be the result of shoaling 
solitons, propagating along the acoustic path. 


Standard Deviation of the In-Situ SIL Estimate (dB) measured over 6 Transmissions (-30 min). 
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Figure 32. The standard deviation of the measured sound intensity (dB) associated with 
the bandpass filtered 10-, 25-, 50-, and 75-Hz bandwidth signal as well as the 
original 100-Hz bandwidth signal. The standard deviation was calculated from the 
measured intensities using a 30 minute window. 
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standard Deviation ofthe In-Situ SIL Estimate (dB) measured over36 Transmissions ("180 min). 
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Figure 33. The standard deviation of the measured sound intensity (dB) associated with 
the bandpass fdtered 10-, 25-, 50-, and 75-Hz bandwidth signal as well as the 
original 100-Hz bandwidth signal. The standard deviation was calculated from the 
measured intensities using a 180 minute window. 


Before proceeding to the next Section, it is important to summarize the results 
from the statistical analysis thus far. First, the statistics of both the modeled and measured 
intensity are clearly sensitive to the bandwidth of the received signal. The variance of the 
sound intensity was observed to increase (decrease) as the bandwidth decreased 
(increased). Second, because the pulse-width of an arrival (i.e., its temporal resolution) is 
inversely proportional to the bandwidth of the received signal, the ability to resolve two 
adjacent arrivals depends not only on their temporal separation, but also on their pulse- 
width (Fig. 31). Thus, based on these two observed results, it can be inferred that as the 
pulse-width of the individual arrivals increases (decreases) due to narrowing (broadening) 
of the received signal’s bandwidth, adjacent arrivals will begin to overlap one another 
(become well separated from one another) resulting in both fewer (more) independent 
arrivals and an increase (decrease) in the magnitude of the variance. These results are in 
complete agreement with the proposed hypothesis. Lastly, the variance associated with 
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the intensity fluctuations are time-dependent, whose magnitude is heavily influenced by 
the periodic presence/absence of a highly energetic internal wave field. However, 
expanding the temporal window of the variance measurements acted to “smooth out” the 
temporal variations of the calculated variance and converged the variance perturbations 
towards a constant value. The convergence of the variance perturbations to a constant 
value suggests that the statistics of the SIL estimates, if calculated over a long enough 
temporal “window,” which would include numerous high frequency wave events, will 
reach saturation. This will be a critical concept in the next Section when validating the 
proposed hypothesis. The final remaining task in Objective 2 is to devise a method that 
will reliably determine the number of arrivals and relate them to the observed variance. 
The methodology used to determine the number of arrivals for each bandwidth will be 
discussed next. 

D. RELATING THE OBSERVED STATISTICS TO NUMBER OF 

ARRIVALS (VALIDATING THE HYPOTHESIS) 

1. Determining the Number of Arrivals 

An objective approach was taken when attempting to determine the number of 
arrivals present in the signal segment. Originally, in order to be consistent with the 
previous Section in this dissertation, the signal segment from 19.0-s to 19.95-s (blue 
segment. Fig. 34) was used when calculating SIL estimates and the noise segment from 
22.0-s to 22.95-s (green segment. Fig. 34) was used for noise analysis. The first task was 
to establish a signal-to-noise threshold that would, when applied, confidently and 
consistently identify arrivals associated with the received signal and not a peak resulting 
from matched-filter noise. This signal-to-noise threshold (green horizontal line. Figs. 34 
and 35) was calculated by finding the peak value in the noise segment, and doubling it. 
This method established a dynamic signal-to-noise threshold based solely on the 
maximum noise level recorded in the noise segment. For a peak to be considered a 
significant arrival, the magnitude of the peak must be +6 dB (i.e., twice the amplitude of 
the peak noise value) above the noise. Arrivals found to exceed this threshold were 
considered significant arrivals originating from the sound source rather than processor 
noise. 
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Time (s) 

Figure 34. Example arrival strueture assoeiated with a 100-Hz bandwidth signal. The 
signal segment used to generate the SIL estimates has been colored blue, while 
the noise segment has been colored green. The green horizontal line is the noise 
threshold which is 6 dB above the highest observed value in the noise segment. 


Examining Eig. 35, one should note the long, temporal width of the first arrival 
(lasting from approximately 19.3-s to 19.42-s). This arrival is obviously composed of 
several overlapping individual ray/mode arrivals and raises the following questions; Do 
the smaller amplitude, later arrivals that exceed the noise threshold contribute to the 
sound intensity level variance? Additionally, what temporal separation should be required 
between two adjacent arrivals in order to be classified as two independent arrivals (i.e., 
“well separated”)? All of these factors were investigated during the course of this 
research and will be discussed in the subsequent paragraphs. 
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Arrival Structure from transmission #3!I0 off of hydrophone 8. 
100 Hz bandwidth. 



Tima (a) 

Figure 35. Same as Fig. 34, however, the signal segment has been enlarged in order to 
show the temporal width of eaeh arrival (red line) within this arrival strueture. 
The peak values within eaeh arrival have also been identified (red asterisk). 


a. Arrivals Found within the Full Signal Segment 

To examine the significance of the late arrivals, two methods were used to 
analyze the signal. The first method considered all arrivals that exceeded the noise 
threshold over the full temporal length of the signal segment (i.e., from 19.0-s to 19.95-s) 
as significant and were used when calculating the DOF. The second method (to be 
discussed in paragraph IV.D.l.b.) would remove many of the later, smaller amplitude 
arrivals by only using the arrivals found to exceed the noise threshold in an abbreviated 
signal segment (i.e., from 19.0-s to 19.4-s) when estimating the DOF. The final question 
to be addressed is assigning the appropriate number of DOF to the arrivals. Although 
other methods were examined (and will be briefly mentioned in the closing paragraph), 
treating broad, overlapping groups of arrivals that were sufficiently separated from 
adjacent arrival groups (i.e., based on the temporal separation requirement) as 

independent degrees of freedom proved the most successful and is discussed next. 
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Using the +6 dB noise threshold as a guide to eonfidently identify 
multipath arrivals from proeessor noise, a method was devised that treated overlapping 
arrivals (henee, arrival groups) as a signifieant, independent arrival. This method 
ealeulated the DOF based on two parameters. The first parameter was simply ealculating 
the number of “arrival groups”. An “arrival group” value of one was assigned to any 
portion of proeessed signal that exeeeded the noise threshold. It simply assigned any 
signal segment that exeeeded the noise threshold two degrees of freedom. For example, 
in Fig. 35 there were six signal segments where the amplitude of the proeessed signal 
exeeeded the noise threshold; eonsequently, twelve degrees of freedom were assigned to 
this arrival strueture. 

The seeond parameter this method attempted to quantify is when adjaeent 
arrivals should be treated as independent arrivals, viee grouped together. For example, in 
Fig. 35, the last two arrivals are very elosely spaeed. If the noise threshold were to 
deerease slightly, then these two arrivals would no longer be separated and the degrees of 
freedom would deerease by two. Consequently, beeause the noise threshold is somewhat 
arbitrary and the temporal position of these arrivals will oseillate slightly about a mean 
travel time over the eourse of several arrival struetures, there eould be ineonsisteney on 
the degrees of freedom assigned to them. To prevent this, a temporal separation 
requirement between adjaeent arrivals was established to aid in the proper DOF 
ealculation. This requirement essentially separated independent arrivals from a group of 
adjaeent arrivals by examining their temporal separation. If the spacing between two 
adjacent arrival groups exceeded the temporal threshold, then the arrival groups were 
treated as two independent arrival groups and were assigned four degrees of freedom. 
Conversely, if the separation between the two adjacent arrival groups did not exceed the 
defined temporal separation, then the two arrival groups were treated as one, broad 
arrival group and were assigned two degrees of freedom. 

Two different definitions of pulse-width were used during this study. The 
first, which will be referred to as the “narrow” pulse-width, was defined as twice the 
temporal width between the peak value and a 3 dB decrease in amplitude. The other, 
which will be referred to as the “broad” pulse-width, was determined by the temporal 
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width of the base of the pulse. Utilizing the two pulse-width definitions as the temporal 
separation requirement, two different thresholds were used in an attempt to isolate 
independent arrival groups. 

The estimated DOF caleulated by the “group” method, utilizing either the 
“narrow” or “broad” pulse-width as the temporal separation requirement, are presented in 
Figs. 36 and 37, respectively. Analysis of Figs. 36 and 37 indicate that the estimated DOF 
experiences severe fluctuations when using either the “broad” and “narrow” temporal 
separation requirement. This was a surprising result and a significant observation of this 
research. The cornerstone of Dyer’s statistical theory and the proposed hypothesis in this 
dissertation is the assumption that the arrival structures are stable and the statistics of the 
intensity fluctuations result primarily from the phase randomness of the received signal. 
Both Figs. 36 and 37 suggest that the arrival structures are not stable, but rather changing 
in time possibly due the high-frequency wave field. In order to determine if the periodic, 
high frequency internal wave field is causing fluctuations in the estimated number of 
arrivals, a comparison between the standard deviation of the estimated DOF and the 
standard deviation of the estimated SIL was made and is discussed next. 

Figure 38 shows the standard deviation of the fluctuations in the SIL 
estimates using a 10 minute window, while Figs. 39 and 40 show the standard deviation 
of the estimated DOF measured over the same temporal range. Examination of Fig. 38 
reveals four broad periods during which the standard deviation of the SIL estimate 
reaches a maximum. Although present in all five bandwidths, the maxima are more 
pronounced in the 50-, 75-, and 100-Hz bandwidth signals. The variance of the SIL 
estimates were observed to increase from 13:00 - 21:00 April 13, 03:00 - 09:00 April 14, 
13:00 - 21:00 April 14, and 05:00 - 08:30 April 14, 2005. 

Comparison between of Fig. 38 and Figs. 39 and 40 reveal a striking 
correlation between high variance periods observed in the SIL estimates to the variability 
observed in the DOF estimates. Regardless of the pulse-width definition, roughly five 
periods of standard deviation maxima in the DOF estimates were observed, namely, 
13:00 - 22:00 April 13, 03:00 - 09:00 April 14, 13:00 - 21:00 April 14, 22:00 - 01:00 

April 14, and 05:00 - 08:30 April 14, respectively. With the exception of the fourth high 
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variance period observed in the DOF estimate time-series, the periods of variance 
maxima between the SIL estimates and DOF estimates almost perfectly coincide. 

Before proceeding to the next Section, it is important to summarize the 
results from the statistical analysis of the DOF estimates and their correlation to the 
sound intensity variance. First, DOF estimates are highly variable. This variability 
indieates that either the arrival struetures are unstable or the noise threshold is fluctuating 
(i.e., the amount of the signal that exceeds the noise threshold is fluctuating which, in 
turn, affects the number of signifieant arrivals (red line. Figs. 34 and 35)). Seeond, 
periods of variance maxima between SIL estimates and DOF estimates appear to 
eorrespond. This is signifieant because the statistics of SIL estimates are not dependent 
upon the noise threshold. Consequently, any observed variability in the intensity 
measurements is completely dependent on the fluctuations of the arrival strueture, not the 
noise threshold used in this Section to detect significant arrivals. Because the standard 
deviation maxima and minima of the DOF estimates and the SIL estimates show sueh a 
strong eorrelation with eaeh other, one can infer that the cyelical high frequeney 
nonlinear internal wave field, rather than the random proeess of ambient or proeessor 
noise, is driving the statistics. Consequently, if noise isn’t eausing the fluetuations 
observed in the DOF estimates, then the shape of the arrival strueture (i.e., both 
amplitude/phase and travel time) must be ehanging. This is a signifieant result. Not only 
are the arrival structures unstable (Figs. 36 and 37), but their variability (i.e., estimated 
DOF) appears connected to the variance induced by the high-frequency nonlinear internal 
wave field. In the next Section, the late arrivals will be removed from the analyzed signal 
and the statistics of the abridged signal are examined. 
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Estimated DOF using the "Group" Method with the "Narrow" Pulse Width Separation Requirement 
(Signal Segment: 19.0-s to 19.9&-s) 

10 Hz Bandwidth Signal 


■I I III II 


Apr.13,2C05 12:CO:00 /Vif-13.2005 18:00:00 Apr.14,2005 00:00:00 Apr.14,2005 06:00:00 Apr.14,2005 12:00:00 Apr 14^05 18:00:00 Apr. 15^005 00:00:00 Apr. 15,2005 06:00:00 

25 Hz Bandwidth Signal 



Apr.13,200S 12:00:00 Apr.13,20IOS 18:00:00 Apr.14,2005 00:00:00 Apr.14,2005 06:00:00 Apt.14,2005 12:00:00 Apt.14305 10:00.00 Apt.152005 00:00:00 Apr.1S,20O5 06-00:00 

50 Hz Bandwidth Signal 



Apr.13,20IO5 12:00:00 Afir.13,20D5 18:00:00 Apr.14,2005 00:00:00 Apf.14,2005 06:00:00 Apr.14,2005 12:00:00 Apr.142005 13:00:00 Apr. 152005 00:00:00 Apr. 15,2005 06:00:00 

75 Hz Bandwidth Signal 



Apr.13,20D6 12:00:00 Apr,13,2005 18:00:00 Apr.14,2005 00:00:00 Apr.14,2005 06:00:00 Apr.142005 12:00:00 Apr. 142005 18:00:00 Apr,152005 00:00:00 Apr.15,2005 06:00:00 

100 Hz Bandwidth Signal 



Apr.13,2005 12:00:00 Apr.13,2005 18:00:00 Apr.14,2005 00:00:00 Apr.14,2005 06:00:00 Apr.14,2005 12:00:00 Apr. 142005 18:00:00 Apr.15,2005 00:00:00 Apr.15,2005 06:00:00 


Time (GMT) 


Figure 36. Estimated DOF using the “group” method and the “narrow” pulse-width 
separation requirement between adjacent arrival groups. 


Estimated DOF using the "Group" Method with the "Broad" Pulse Width Separation Requirement 
(Signal Segment: 19-sto 19.95-s) 

10 Hz Bandwidth Signal 



Apr,132005 12:00:00 Apr.13,2005 18:00:00 Apr 143005 00:00:00 Apr.14,2005 06:00.00 Apr.14 2005 12:00:00 Apr.1420D5 18:00:00 Apr.15 2005 00:00:00 Apr,152005 06:00:00 

25 Hz Bandwidth Signal 
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£ 
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Apr.133005 12:00:00 Apr.13,2005 18:00:00 Apr.142005 00:00:00 Apr.14,2005 06:00:00 Apr.14,2005 12:00:00 Apr. 14,2005 18:00:00 Apr.15,2005 00:00:CD Apr. 152005 06:00:00 

50 Hz Bandwidth Signal 



Apr.132005 12:00:00 Apr.13,2005 18:00 00 Apr 14 3005 00:00:00 Apr.14,2005 06:00.00 Apr 14 3005 12:00:00 Apr.142005 18:00:00 Apr 15 3005 00:00:00 Apr 15 3005 06:00:00 

75 Hz Bandwidth Signal 



Apr.13 2005 12:00:00 A^r 13,2005 18:00:00 Apr.142005 00:00:00 Apr.14,2C05 06:00:00 Apr.14,2005 12:00:00 Apr.14,2006 16:00:00 Apr.15,2005 00:00:00 Apr.152005 06:00:00 

100 Hz Bandwidth Signal 



Apr.132005 12:00:00 Apr. 13,2005 18:00:00 Apr.142005 00:00:00 Apr.14,2005 06:00:00 Apr.14,2005 12:00:00 Apr.14,2005 18:00:00 Apr.16,2005 00:00:00 Apr.152005 06:00:00 

Time (GMT) 

Figure 37. Estimated DOF using the “group” method and the “broad” pulse-width 
separation requirement. Note the range scale is different from that of Fig. 36. 
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Depth (m) Depth (m) Depth (m) Depth (m) Depth (m) 


Standard Deviation of the In-Situ SIL Estimate (dB) measured over2 Transmissions (-10 min). 

10 Hz Bandwidth Signal 

Illll.lll nil mil lll■llll BM 





Apr.13,2005 12;00;00 Apr.13,2005 18:00:00 Apr.14^006 00:00:00 Apr.142005 06:00:00 Apr.142005 12:00:00 Apr.U,2005 10:00:00 /^)r,15,2005 00:00:00 Apr.15^005 06:00:00 

25 Hz Bandwidth Signal 

miniiuii 




..iJMum'Li.r.imi.. 41.11 imujiii.... JiJU' 'i.' 




Apr.13,2005 12:00:00 Apr. 13,2005 18:00:00 Apr 14,2005 00:00:00 Apr. 142005 06:00:00 Apr.142005 12:00:00 Apr.14,2005 18:00:00 Apr.15,2005 00:00:00 Apr. 15,20[S 06:00:00 

50 Hz Bandwidth Signal 





Apr.13,2C05 12:00 00 Apr.13,2005 18:00:00 Apr.14,2006 00:00:00 Apr. 142005 06:00:00 Apr.142005 12:00:00 Apr.U,20CK 18:00:00 Apr.15,2005 00:00:00 Apr.15,2006 06 00:00 

75 Hz Bandwidth Signal 


IIL ..ldkAiUiJ'4U!L;.i j;.l!:llllU4iil«.J.JLVi 



Apr.13,20D5 12 00:00 Apr.13,2005 18:00:00 Apr.14 2005 00:00 00 Apr 142005 06:00:00 Apr.1420105 12:00:00 Apr.14,2005 18:00:00 Apr.15,2005 00:00:00 Apr.1520C6 06:00 00 

100 Hz Bandwidth Signal 

lllini'iil lTTAWiriHBll.JjaiHI'iUL .il4IIUUIIIilJ'.IIIIUiJ.i 4111 
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: 1 m 
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Apr.13,2005 12:00:00 Apr.13 2005 18:00:00 Apr,14,2006 00:00:00 Apr.U^OOS 06:00:00 Apr.142005 12:00:00 Apr,142005 18:00:00 Apr,15,2005 00:00:00 Apr.l5 2005 06:00:00 

Time (GMT) 

Figure 38. Standard deviation associated with the measured SIL estimates (dB) as 
measured over 10 minutes. 
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Standard Deviation of the Estimated DOF using the "Group" Method with the "Narrow" Pulse Width Definition 
Measured over 2 Transmissions (~10 minutes). 

10 Hz Bandvfldtii Signal 


■ I I I I IIIH II II 



Apr.13.2005 12:00:00 Apr.13.2005 18:00:00 Apr.142005 00.00:00 Apr.142005 06:00:00 Apr.14,200512:00:00 Apr.14.2005 18:00:00 Apr.15.2005 00:00:00 Apr.15.2005 06:00:00 

25 Hz Bandwidth Signal 


I 



Apr.13200S 12:00:00 Apr.132005 18:00:00 Apr,142005 00:00:00 Apr.142005 06:00:00 Apr.14,2005 12:00:00 Apr 14,2005 18:00:00 Apr,15,2005 00:00:00 Apr.152005 06:00:00 

50 Hz Bandvfldth Signal 



Apr.132005 12:00:00 Apr.132005 18:00:00 Apr.142005 00.00:00 Apr.142005 06:00:00 Apr. 14,2005 12:00:00 Apr.14,2005 18:00:00 Apr.15,2005 00:00:00 Apr.152005 06.00:00 

75 Hz Bandwidth Signal 

.t.r U l.l I 



Apr.13,2005 12:00:00 Apr.13200518:00:03 Apr.142005 00:00:00 Apr.14,2005 06:00:00 Apr.U,2005 12:00:00 Apr.14,200510:00:00 Apr.15,2005 00;CD:00 Apr.lS^OOS 06:00:00 

100 Hz Bandwidth Signal 

iitmUllfllliiiiiiU Am r:ii |: ii I 

Apr.13,2005 12:00:X Apr.132005 18:00:00 Apr.142005 00:ra:00 Apr.14,2005 06:00:00 Apr.14,200512:00:00 Apr.14,2005 18:00:00 Apr.15,2005 00:00:00 Apr.152005 06:00:00 

Time (GMT) 

Figure 39. Standard deviation of the estimated DOF from Fig. 36 (i.e., “group” method 
and “narrow” pulse-width). 



Standard Deviation of the Estimated DOF using the "Group" Method with the "Broad" Pulse Width Definition 
Measured over 2 Transmissions (^10 minutes). 

10 Hz Bandwidth Signal 



Apr.13,2005 12:00:00 Apr.13,2005 18:00:00 Apr.14,20C6 00:00:00 Apr.14,2C05 06:00:00 Apr. 14,3305 12:00:00 Apr.14,2005 18:00:00 Apr.153305 00:00:00 Apr.15,3305 06:00 OO 

25 Hz Bandwidth Signal 



Apr.13,20D5 12:01:00 Apr.132005 18:00:00 Apr.14,20D5 00:00:00 Apr.14,2005 06:00:00 Apr. 143005 12:00:00 Apr.14,2005 18:00:00 Apr. 152005 00:00:00 Apr. 152005 06:00:00 

50 Hz Bandwidth Signal 

T 


"i'll' fl' 

; mil i ni 



1 ! I 




Apr.13,200512:00:00 Apr.13,2005 10:00:00 Apr.14,2005 00:00:00 Apr.14,2005 06:00:00 Apr.14,2005 12:00:00 Apr.14,200518:00:00 Apr.152005 00:00:00 Apr.1523D5 06:00:00 

75 Hz Bandwidth Signal 

nifiriH'iiii" rm, v "t 





Apr.13,2005 12:00:00 Apr.13,2005 18:00:00 Apr.14,2005 00:00:00 Apr,14,2005 06:00:00 Apr.14,2005 12:00:00 Apr.14,200518:CD:CD Apr.152005 00:00:00 Apr.16,3305 06:00:00 

100 Hz Bandv^dth Signal 

II ■Hill rniH Tu 

Apr.13,200S 12:00:00 yM3r.13,2005 10:00:00 Apr.1420C600:0000 Apr.14,20ffi0e:00:00 Apr.14200512:0000 Apr.14,200S 18:00:00 Apr.15 2005 00:00:00 Apr 15200506:00:00 

Time (GMT) 

Figure 40. Standard deviation of the estimated DOF from Fig. 37 (i.e., “group” method 
and “broad” pulse-width). 
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b. Estimating the Number of Arrivals Using other Methods 


As mentioned before, an “abridged” portion of the signal was examined in 
order to analyze the significanee of the late arrivals. The “abridged” signal was defined 
from 19.0-s to 19.4-s. The 19.4-s cutoff was selected because it was found to be a good 
separation point between the early, high amplitude arrivals and the later, small amplitude 
arrivals. It should be noted that the cutoff of 19.4-s was not fixed. In some cases, arrivals, 
whose amplitude exceeded the noise threshold, would extend beyond the 19.4-s cutoff 
As long as more than 50% of its arrival width was found before the 19.4-s cutoff, the 
arrival was used to estimate the DOF. As before, both the “narrow” and “broad” pulse- 
width definitions were used as the temporal separation requirements between adjacent 
arrival groups and the results are discussed next. 

Similar to results in the previous Section, the estimated degrees of 
freedom using the abridged signal segment fluctuated with time. Variability maxima in 
the DOF continued to coincide with variance maxima from the intensity estimates (not 
shown), again suggesting that the high-frequency nonlinear internal wave field was 
dominating the short-term statistics. Due to the “abridged” signal and temporal separation 
requirement between arrivals, a minimal number of DOF were calculated using both the 
“narrow” and “broad” temporal separation requirements. Because the main arrival group 
was very broad and dominated the “abridged” signal, this method significantly 
underestimated the estimated DOF. 

2. Relating Average Number of Arrivals to the Measured Variance of 
the Sound Intensity Level 

There are three critical points to understand how the proposed hypothesis relates 
the number of arrivals to observed variance of the intensity measurements. First, this 
hypothesis assumes that the arrival structure is stable (i.e., the number of non-overlapping 
arrivals is invariant or, if the estimated number of arrivals does fluctuate, then the 
deviation of the number needs to be small). Second, the phase randomness of each 
individual arrival resulting in either constructive or destructive interference of the 
pressure field when reaching the hydrophone concurrently with other arrivals will drive 
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the observed statisties of the intensity measurements. Third, the hypothesis requires a 
varianee time window large enough to colleet a suffieient number of samples from 
numerous wave events in order to reach phase randomness in the sample population. 

The results from the estimated DOF study clearly suggest that arrival structures 
were not stable along the transmission path. The number of arrivals fluctuated as a result 
of the passage of nonlinear internal waves along the acoustic path. As a result, for each 
variance measurement of the acoustic intensity, the estimated number of arrivals 
fluctuated. Consequently, within each of these variance windows, there existed a sample 
population of estimated number of arrivals. Assuming the proposed hypothesis is valid 
and the three critical points are met, this sample population should be able to relate the 
estimated number of arrivals to the measured variance of the acoustic intensity (i.e., 
follows the curve in Fig. 28). Note that for each variance measurement of the sound 
intensity, there is a corresponding estimated DOF sample population distribution where 
both the mean number of arrivals and the variance of the number of estimated arrivals are 
always greater than zero. The fact that the estimated DOF has a distribution is an 
important concept. As with any distribution, there will be the most common value (i.e., 
the mean number of estimated DOF) and other moments that can be calculated. Thus, for 
each variance measurement of the sound intensity, there will be corresponding mean 
number of arrivals determined. Attempts to relate the most common (i.e., averaged 
estimated DOF) value to the observed intensity fluctuations will be discussed next. 

Figure 41 depicts the evolution of a series of scatter plots based on the average 
number of arrivals and the corresponding variance in the acoustic intensity using 
progressively larger variance windows for four out of the five broadband signals. In other 
words, for each variance measurement of the acoustic intensity, an average estimated 
DOF (i.e., twice the averaged number of arrivals) was calculated from the sample 
population collected over the defined variance window. Once the averaged estimated 
number of arrivals and variance of the observed intensity was calculated, they were 
plotted against each other in Fig. 41. It should be noted that when the results from the 75- 
Hz and 100-Hz bandwidth signal were plotted, they tended to overlap each other. In this 
study, the received bandwidth of the 400-Hz sound source signal was assumed to be 100- 
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Hz. However, the true bandwidth of the reeeived signal was probably eloser to 75-Hz. 
This would explain why the results between the two broadband signals overlapped one 
another. Because plotting both the results from both broadband signals does not add any 
information to the analysis, only the results from one of the signals will be presented in 
the scatter plots (i.e., the original 100-Hz bandwidth signal). 

If the proposed hypothesis is valid, the method of estimating the number of 
arrivals is correct, and the three critical requirements are met, then the scattered values 
from each of the broadband signals should converge toward the red line depicted in Fig. 
28. Using the methodology from the “full signal” study and a temporal separation 
requirement equivalent to 80% of the “narrow” pulse-width definition, the average 
number of resolvable arrivals and the associated standard deviation of the acoustic 
intensity for four expanding variance windows were calculated and plotted. Starting in 
the upper left and working clockwise. Fig. 41 depicts how step increases in the variance 
window affects the distribution of the statistics of both the measured intensity variance 
and the average estimated number of arrivals. By increasing the variance window from 
10 minutes to 90, 300 and finally 1440 minutes (i.e., 24 hours), the variance of the sound 
intensity for each broadband signal converged toward a constant value (i.e., less 
scattering along the vertical axis). Additionally, the larger sample size of the estimated 
number of arrivals drove the mean of the sample population toward the true mean of the 
parent population (i.e., less scattering along the horizontal axis). The agreement between 
the observed variance of intensity and the average number of arrivals found within the 
24-hour variance window is excellent. 
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Figure 41. When the variance time window is extended temporally, the proposed 
hypothesis that relates the average number of arrivals to the statistics of the 
acoustic intensity converges toward the theoretical curve. Similar to Fig. 24, the 
horizontal axis is the averaged estimated number of arrivals while the standard 
deviation associated with the measured SIL (dB) is plotted on the vertical axis. 
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Taking into consideration the varianee windows ehosen when ealculating the 
standard deviation of the intensity fluctuations and the mean of the estimated number of 
arrivals, the results analyzed in Fig. 41 were expeeted. When using a short-term varianee 
window (i.e., 10 and 90 minutes), the scatter plots depiet the attempt to relate the average 
number of arrivals to a statistieal measurement that does not meet the three eritical 
assumptions of the proposed hypothesis. As observed in the previous Seetion, when using 
sueh a short time window to estimate the DOF, the arrival struetures were not stable and 
the varianee of the estimated number of arrivals was large due to the presenee/absenee of 
highly energetic, nonlinear internal waves along the acoustic path. These temporally 
variant internal wave aetivities add significant variability to the short-term statisties of the 
acoustic intensity measurements and results in the vertical spread of their measured 
statisties in the seatter plots. Also consider that, at most, a 10 minute varianee window 
only provides 448 independent realizations of the estimated DOF (assuming all eight 
hydrophones didn’t experience any hardware failure). With sueh a limited sample size, 
the sample population will not eontain enough samples from nonlinear wave events to 
ensure that the sample statisties originate from a sample population whose phase 
distribution is uniformly random. With only 448 realizations, the sample population of 
the estimated DOF was too small to eonfidently determine the mean of the parent 
population (i.e., the true DOF) as well. By expanding the variance window, the sample 
size will increase, resulting in a sample population whose statistieal distributions will be 
based on a uniform phase random distribution and whose number of temporally disjoint 
arrivals will approaeh a constant ensemble mean while the variance of the estimated 
number of arrivals approaehes zero. 

One last examination of the Fig. 41 reveals a series of roughly four “varianee 
bands” when the varianee window was expanded to 1440 minutes (i.e., 24-hours). This 
banding can be observed in the seatter plot of all four bandwidths. These bands resulted 
from the teehnique used to ealeulate the standard deviation of the measured SIL. Because 
all eight hydrophones were used when caleulating the varianee, the sample population 
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from a hydrophone that experience intermittent failure (i.e., hydrophones 1 and 3) had 
slightly different statistics than a hydrophone that did not experience hardware failure 
(i.e. the bottom three hydrophones). 
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V. CONCLUSIONS 


The sound-speed field found along the shelf region in the northeast portion of the 
South China Sea is extremely dynamie where varying temporal and spatial seale 
proeesses dominate the observed fiuetuation in the sound-speed field. Signifieant 
variability in the sound-speed field was observed to be indueed by the broad 
expansion/eompression of sound-ehannel assoeiated with the nonlinear internal tides as 
well as the narrow, high-frequeney, nonlinear internal depression and elevation waves 
superimposed on this long-wave pattern. 

Analysis of the measured aeoustie transmissions revealed that the mean vertieal 
distribution of the sound intensity time-series followed the eompression/expansion of the 
sound-ehannel assoeiated with the internal tides. Additionally, the strong horizontal 
sound-speed gradient found along the leading (trailing) edge of elevation (depression) 
waves effeetively seattered sound energy into higher (higher and lower) aeoustie modes. 

Two different vertieal distributions of the sound intensity fluetuations were 
observed depending on the type of episode. Standard deviation maxima of the measured 
SIL were eonfined to the lower half of the water eolumn when elevation wave paekets 
were observed in the vieinity of the hydrophone array. With a eompaeted sound-ehannel, 
the seattering of aeoustie energy was limited to very low aeoustie modes. Additionally, 
energy that was seattered into higher modes was stripped off after several interaetions 
with the bottom, thus making it an extremely loealized event. Additionally, when 
depression waves were transiting over the sound souree, the variability in the SIL was 
eonfined to the lower half of the water column due to the suppression of higher acoustic 
modes near the sound source. With most of the sound energy initially limited to the 
lowest acoustic modes, acoustic energy could only be scattered out of the lowest acoustic 
modes. However, when depression waves transited across the hydrophone array, the 
standard deviation maxima of the observed SIL was much more vertically diffuse. During 
these events, the expanded sound-channel driven by the internal tides allowed for the 
excitation and propagation of additional acoustic modes. With a larger number of 
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acoustic modes available to seatter energy, the strong horizontal sound-speed gradient 
found along the trailing edge of these narrow, depression waves effeetively seattered 
aeoustie energy throughout the water eolumn. 

EOF deeomposition of the temperature time-series eolleeted from the five 
environmental moorings allowed for the generation of an empirieal, two-dimensional 
sound-speed field that was eontinuous in both spaee and time for use within an aeoustie 
model. Comparisons between the modeling results from the eoupled, normal-mode 
aeoustie model and the analyzed data from the hydrophone array revealed that the 
empirieal sound-speed model was able to reasonably reproduee the vertieal distribution 
of the sound intensity time-series resulting from the long-wave pattern, as well as the 
vertieal distribution of the sound intensity fluetuations resulting from pulse propagation 
through an aeoustie path oeeupied by nonlinear depression waves. The empirieal sound- 
speed model, however, was not able to replieate the transformation of a depression wave 
into a series of elevation waves. Consequently, beeause the empirieal sound-speed model 
eould not model the transforming nonlinear internal waves observed in the mooring data, 
no explanations toward the phenomenology of the observed sound intensity fluetuations 
resulting from elevation waves near the reeeiver ean be made. However, the eonfinement 
of the vertieal distribution of aeoustie energy fluetuations to the lower portion of the 
water eolumn may have been due to the eompressed sound-ehannel assoeiated with an 
elevation wave event. With a eompressed sound-ehannel, only the lowest aeoustie modes 
would be able to propagate along the aeoustie path, resulting in most of the aeoustie 
energy foeused toward the bottom. The strong horizontal sound-speed gradient assoeiated 
with the edge of the nonlinear elevation waves may have allowed weak seattering of the 
bottom-foeused, aeoustie energy into higher neighboring aeoustie modes; however, 
beeause of the eompressed sound-ehannel, the weak aeoustie energy seattered into these 
higher aeoustie modes may have been quiekly been stripped away by the bottom. 

The aeoustie modeling results during whieh time a depression wave was observed 
along the aeoustie path helped to explain some of the measurements noted from the 
hydrophone data. First, the eompression/expansion of the sound-ehannel along the 
aeoustie path (i.e., the long-wave pattern) eontrolled both the vertieal distribution of the 
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acoustic intensity and which acoustic modes were initially excited. Although aeoustie 
modes from three acoustic mode groups were initially excited when in elose proximity to 
the sound souree, only the lowest members (1-3) of the low aeoustie mode group were 
able to reaeh the hydrophone array when the sound-ehannel was compressed, as higher 
aeoustie modes were either suppressed or quiekly stripped away by the bottom. As the 
sound-channel expanded over either the partial or the full length of the aeoustie path, 
additional members (1-6) of the low aeoustie mode group were initially exeited and 
aeoustie energy spread into the mid- and lower portions of the water eolumn. Seeond, 
severe fluetuations in the vertical distribution of the sound intensity appear to result from 
the strong horizontal sound-speed gradient found along the leading (trailing) edge of 
elevation (depression) waves. The vertical scattering of sound energy was linked to 
strong aeoustie mode-coupling along the nonlinear internal wave-front. The acoustic 
modeling results suggest that the nonlinear internal wave-front interfaee allows for the 
effeetive seattering of acoustic energy from low aeoustie modes into higher aeoustie 
modes. However, if the scattering of energy into higher aeoustie modes due to acoustic 
mode-coupling occurs too far from the hydrophone array, the aeoustie energy scattered 
into the higher aeoustie modes may be stripped out by the bottom. Finally, although three 
aeoustie mode groups are initially excited near the aeoustie source, the intermediate and 
high aeoustie mode groups appear to eontribute little (if any) energy to the sound 
intensity level. Acoustic energy from the highest aeoustie mode group (18-25) is quickly 
stripped by the bottom prior to reaehing the hydrophone array, while the intermediate 
aeoustie mode group (8-16) appears to eontribute minimally to the observed sound 
intensity. 

The results from the seeond objeetive indieate that is possible to prediet the 
varianee of the transmission loss based on the mean number of arrivals if the temporal 
separation between adjaeent arrivals is properly defined and the phase distribution of the 
sample population is uniformly random. As indieated in the hypothesis, the number of 
resolvable arrivals was found to depend on the reeeived signal bandwidth and the criteria 
of “well-separateness.” Once a teehnique was established to properly identify significant 
arrival against processor noise, four different methods of ealeulating the number of 
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arrivals were examined to assess the significanee of early versus late arrivals (i.e., full vs. 
abridged signal length), as well as two different temporal separation requirements 
between adjacent arrival groups based on two definitions of pulse-width. The estimated 
number of arrivals was observed to fluctuate wildly in all four methods, thus indicating 
an unstable arrival structure. A strong correlation was observed between the standard 
deviation maxima of the estimated number of arrivals and the standard deviation maxima 
of the measured sound intensity suggesting that the high frequency nonlinear internal 
wave field, not the background noise, is driving the variability observed in the short-term 
statistics. Attempts at relating the estimated number of arrivals to the short-term 
fluctuations in the sound intensities were unsuccessful for several reasons. First, the 
proposed hypothesis requires that the standard deviation of the estimated number of 
arrivals be small. When using a short variance window, the standard deviation of the 
estimated number of arrivals will be large, especially when highly energetic, nonlinear 
internal waves are found along the acoustic path. If the variance window is expanded, the 
sample size will increase, thus reducing the variance of the estimated number of arrivals. 
Additionally, as the sample population is increased, the confidence in the mean of 
estimated number of arrivals will increase as the sample mean approaches the parent 
mean. Second, the proposed hypothesis requires a sample population whose statistical 
distributions are based on a uniform phase random distribution. Short-term statistics 
provide a limited sample size, thus the sample population will not contain enough 
samples from the various nonlinear wave events to ensure that the calculated statistics 
originate from a sample population whose phase distribution is uniformly random. If the 
variance window is expanded, the sample population increases and this discrepancy can 
be overcome. Lastly, the short-term standard deviation of the measured sound intensity is 
time-variant due to the presence/absence of highly energetic, nonlinear internal waves 
along the acoustic path. These temporally variant internal wave activities resulted in a 
significant increase in the volatility of short-term statistics of the acoustic intensity 
measurements. Again, by expanding the variance window, the temporal variability in the 
standard deviation of the observed sound intensity was minimized. 
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Once the varianee window was allowed to expand to a neeessary temporal width 
(i.e., approximately 12- to 24-hours) and the temporal separation requirement between 
adjacent arrival groups was properly defined, the alignment between the theoretieal eurve 
of the proposed hypothesis and the seatter plots which related the estimated number of 
arrivals to the standard deviation of the sound intensity was excellent. 
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